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SUMMARY 


This report is the first supplement to the primary HILTOP program 
document published in January, 1873, and describes the modifications and im- 
provements made to the HILTOP electric propulsion trajectory optimization 
computer program up through the end of 1974. 

New program features include the simulation of power degradation, 
housekeeping power, launch asymptote declination optimization, and powered 
and unpowered ballistic multiple swingby missions with an optional deep space 
burn. 

The report contains the new analysis describing these features, a com- 
plete description of program input quantities, and sample cases of computer 
output illustrating the new program capabilities. 
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I. 


INTRODUCTION 


This document describes the modifications and improvements made to 
the HILTOP electric propulsion trajectory optimization computer program [l 3 
during calendar years 1973 and 1974 c New program features include the simu- 
lation of power degradation, housekeeping power, launch asymptote declination 
optimization, and powered and unpowered ballistic multiple swingby missions 
with an optional deep space burn. 

The power degradation model has been hypothesized by the authors in 
earlier publications [2, 3l . The model allows a single parameter (denoted 
"characteristic degradation time") to describe the power degradation behavior 
of an electric propulsion spacecraft to a degree which fundamentally affects 
the solution to the trajectory optimization problem. 

The option of simulating spacecraft housekeeping power applies to solar 
electric propulsion with specified reference power. The housekeeping power 
is a specified constant power generated by the solar arrays and shunted away 
from the thruster power-conditioners and directly to the spacecraft payload 
for "housekeeping" purposes. 

The launch asymptote declination optimization model was first hypothesized 
by the authors in the appendix of [lJ , and later a more thorough treatment of the 
subject was put forth in [4] . A solution to the problem of optimizing electric 
propulsion heliocentric trajectories, including the effects of geocentric launch 
asymptote declination on launch vehicle performance capability, is developed 
using variational calculus techniques. The model of the launch vehicle per- 
formance includes a penalty associated with a non-easterly launch plus another 
penalty arising from a non-coplanar launch from the parking orbit. Provisions 
for range safety constraints are included. Optimal trajectories will generally 
have the launch excess velocity offset from the initial primer vector. 
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Extra-ecliptic mission simulations now involve launches from the Earth 
in which the Earth's ephemeris is generated by the program’s analytic ephemeris 
capabilit}'; previously, extra-ecliptic missions were generated simply by starting 
the trajectory on the x-axis at one astronomical unit from the sun. This improved 
extra-ecliptic capability allows the launch date to be optimized together with the 
launch asymptote declination. Also, an additional set of boundary conditions has 
been added to the program for simulating extra-ecliptic missions. 

The program's ability to simulate all-ballistic missions has been ex- 
panded to include powered and unpowered multiple swingby missions with an 
optional deep space burn. This capability renders HILTOP a powerful tool for 
ballistic mission design and optimization, with tremendous flexibility for creating 
imaginative multi-target mission profiles. 

This report contains the analysis describing these new program features, 
a complete description of program input quantities, and sample cases of computer 
output illustrating the new program capabilities. 
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II. 


FORMULATION 


A. SPACECRAFT AND TRAJECTORY MODELS 


The assumed spacecraft and trajectory models are as described in Cl J 
and are not repeated here; the nomenclature used in the following analysis is 
also described in [l] , except for the introduction of new symbols which are 
described in the text as they appear. 


1. Power Degradation . Historically, this electric propulsion power 
degradation model first appeared in the literature in [2] , and then soon after- 
ward an improved discussion appeared in [3] , in which several of the ramifi- 
cations and consequences of the theory are also discussed. For deeper insight 
into the analysis, the reader is therefore referred to [3] , from which much of 
the anatysis below is extracted. 

The model discussed here is concerned with the manner in which the 
performance of a solar array degrades due to high energy particle damage. 

It is assumed that one can define a damage factor q which has a value of unity 
at the initial time and which decreases in value with time during the course of 
the mission such that the power output, p, of the arrays at any time may be 
written 


P 


= q yp 


ref 


where p ^ is the reference power and y is the power factor which is a 
function of solar distance and array orientation relative to the sun. The damage 
factor q may also be thought of as a time -dependent efficiency. The derivative 

o 

q is negative, and is assumed to be linearly proportional to q and to the density 
of particles impinging on the solar cells. For simplicity, it is assumed that the 
particle density d is of the same form as the density of photons striking the 
surface of the array, i.e. , 

e • n 

d= -V « 

r 
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where e^ is a unit vector along the sun-spaeecraft line, n is a unit vector 

normal to the a r rays such that e * n > 0 implies the front of the panels faces 

r 

the sun, and r is the solar distance of the spacecraft,, In the discussion to 
follow, an upper case symbol will denote a vector, a lower case symbol will 
denote a scalar, and a lower case symbol with a bar will denote a unit vector. 



q = - kqd 

(2) 

and 

t 



-k f d dt f 



q = e J Q 

(3) 


where k is the constant of proportionality. It is convenient to introduce a 
parameter s, called degradation time. 


or 


t 

s = J ddt* 
0 

c 

s = d 


(4) 

(5) 


Thus, under the above assumptions the degradation time accumulates at a faster 

rate when the spacecraft is nearer the sun, which is a characteristic one might 

expect. Note that when the spacecraft is maintained at 1 AU with the panels 

normal to the sun line, s = 1 and the degradation time is equal to the flight 

time. Then k may conveniently be thought of as the inverse of a reference time, 

called the characteristic degradation time denoted by T , and the degradation 

d 

factor q becomes 


-s/T. 
q = e d 


(6) 


Actually, there is little reason to allow the degradation time s to continue to 
accumulate during coast phases since the arrays may be turned edgewise to the 
sun and the degradation process may be effectively halted. Therefore, we ad- 
just equation (5) to read 

s = h d (7) 

cr ' ' 
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where h^_ is a step function equal, to one during thrust phases and equal to zero 
during coast phases. 

The characteristic degradation time T is an engineering parameter that 
must be determined experimentally. For example, by exposing a solar cell to 
the particle emission of a solar simulator and measuring the performance of the 
cell over a period of time, one should be able to estimate a reasonable value of 
T . Another source of information would be measurements from actual space- 
craft which employ solar cells for power supply. 

The assumed exponential form of the degradation factor, although intended 
for use with SEP systems, i.s applicable for NEP systems as well. The principal 
difference is in the definition of s ; for example, s’ = h . The exponential form 
also permits one to evaluate radio-isotope systems by defining s = 1 and letting 
T d represent the time for the radioactivity to dissipate to l/e of its initial 
level . 

In the development which follows, the formulation applicable to SEP is 
used exclusively. The equations of motion and adjoint equations are given in 
[l] ; the equations affected by the inclusion of power degradation in the model 
are given below. 

Power degradation affects the problem in a very fundamental sense, be- 
ginning with the rocket-thrust term in the equations of motion: 

V = K = a V 

r (8) 

R = V 

where R is the position vector, r is the magnitude of R, V is the velocity 
vector, /i is the gravitational constant of the sun, a is the magnitude of the 
thrust acceleration and e ^ is a unit vector in the direction of thrust. The 
thrust acceleration a is a function of several variables and may be written as 
follows : 
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a = h 

cr 


gyq 

v 


where g is a reference thrust acceleration evaluated under a prescribed set of 
conditions, v is the ratio of current to initial mass, and q is the degradation 
factor defined above. The power factor y is assumed to be a function of the 
density, d, of photons incident on the arrays, where d is as written in equation 
(1). The general form of y is 
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y = d ) 
i=0 



(9) 


The coefficients a. are chosen so that this equation will adequately describe the 
power output of a given array. The only restriction on the a. is that their sum 
is equal to one. Then at r = 1 AU, with the arrays normal to the sun line, 
y = d = 1. The summation term in (9) represents the temperature effect. 

The mass ratio satisfies the differential equation 

; = -h sjn., do) 

cr c 7 


using v = 1 as an initial condition, where c is the jet exhaust speed which is 
assumed to be constant over the trajectory. 

The variational Hamiltonian becomes 

h = A° R—A • R + A v + A f+As’+Ag+Ac 
v v T s g c 

= h f-^SL (A • i - - A )+A d+A 1 - ^ (A • E) - A • E (11) 
a L p v t c v' s t J 3' ’ ' ' 

r 


and the adjoint equations are 


a=_4a + ^(r. A)R+h (A- 5 - v -\ )+ x 1 

3 5 ' ' aL v t c v 7 sJ 


ad 

br 


A* = h 


gyq 


(A • e.) 
v a .2 ' t 7 


( 12 ) 
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X 

g 


yq 

i — 

a v 


(A * e 


“ * ) 
c v 


where 


A* 

c 



A 

v 



.m (A .- _z x , 

v t e y 


4 

7* = By/dd = 2, ^ ( 1 + ^ 

i=0 


dd 

8R 



n - 3(e 

r 



(12) 

(cont 


The control variables are the thrust direction e , the switch step 
function h , and, providing the array orientation is not constrained to yield 
maximum power, the normal direction n. According to the Maximum Principle, 
these controls are chosen to maximize the variational Hamiltonian (11). The 
maximum of li^ with respect to h^ is seen to depend totally on the sign of the 
term in square brackets. That is, denoting 


a _ JL22L (A e 

V 


V 


, - - A ) + A d + A 
t c v’ s T 


(13) 


then choose 


_ / 1 if CT > 0 

a to if ct < 0 ' 


(14) 


Maximizing h^ with respect to e^ is also accomplished by inspection. Since 
e^ appears only in the dot product with A and since the coefficient of that dot 
product is non-negative, i.e„, 


h^gyqA' s 0 


then h^ is maximized with respect to e^. by making A * e^ as large as possible, 
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which is the same result as when there is no degradation* The control vector 
n appears explicitly in (11) through the density d as given bj' equation (1). In 
fact, siiice n only appears in h through d, n only affects through its dot 
product with e_^. Letting the angle between n and e^ be denoted X such that 

e • n = cos X > then it is clear, that there will be a "best" angle X between n 

r 

and e^ to maximize h , but that n may lie along any element of a right circu- 
lar cone of half angle X about e . For the moment, we will put aside the question 
of the explicit direction of n and concentrate on defining the optimum x or, 
alternatively, the optimum d. The optimum value of d is determined by maxi- 
mizing b^ with respect to d, i.e. , by solving for the root of the equation 

Sr = h a!f ' 0 (“> 


Performing the indicated differentiation yields 


a + X = o 

r s 


y* — - 


where, using (9), 


4 

/* = Y 


V ... i . ,i/4 

2 V 1 + i )d 

i=0 


ct = A • e - - Pi 

r t c V 

i/4 

Because of the form of (18) equation (16) is a quartie in the variable d , and 
is solved by iteration in the program. A more detailed discussion of the solution 
to (16) is given in reference [3] for a specific set of coefficients a^. For now, 
assume that the optimum value of d is found from (16). Then the optimum angle 


X is immediately obtained 


8 


cos X = dr 


2 


(19) 


Of course, equation (16) does not take into consideration the fact that d can 

never exceed the inverse square of the solar distance. Consequently, the right 

hand side of (19) may exceed unity, under which condition the program sets 

cos X 5 1 (i.e., X“0) impljung that n is directed along e . 

r 

If the d that represents the solution to (16) exceeds the upper limit 

allowed for d whether that limit is imposed by the problem or by nature (i.e. , 

/ 2 

1/r ), the correct choice for d is that upper-limiting value. Likewise, on the 
lower side, d is physically limited to be non-negative. Therefore, a negative 
solution to (16) is disregarded, and d is set to zero which corresponds to 
X - tt/ 2 (panels oriented edgewise to the sun), and the engines are shut down. 

The precise definition of n has no bearing on the solution of the problem, 
except as it affects d as defined in (1). The appearance of n in the state and 
adjoint equations and the variational Hamiltonian is solely through the density d 
except in the equation for A where n appears explicitly as part of the partial 
dd/dR. defined following equations (12). Actually, this partial is valid only if 
d is permitted to vary with R. That is, if either d = 0 or d = constant is im- 
posed then dd/3R becomes the null vector and the entire term drops from the 
equation for A . Furthermore, if d is the solution to (16), then the last term 
of A in (12) again drops out because the term in square brackets is the left side 

of (16). Therefore, the only time the term in question remains in the equation 

.. 2 — - 
for A is when d = l/r which corresponds to cos X = e r • n = 1 and implies 

n = e . Under this condition 
r 

dd _ 2_ - 

dR " 3 e r * 

r 

The boundary condition pertaining to the initial degradation time is 


s (t ) = 0 
o 
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Although transversality conditions are discussed in the next section, the trans- 
versality condition associated with degradation time is stated here for the sake 
of completeness. If the final degradation time, s^, is unspecified, 

X = 0 . 

S f 

The initial value of X is unknown and therefore becomes one of the independent 

s 

parameters of the boundary value problem. From equations (12), it follows that 


o 

X 

s 



and, therefore, 



(X - X ) 

' rr 


g. 


g. 


Using the boundary conditions X = X = 0, it follows that 

s f S Q 



Thus, if an approximate value of X is available from a trajectory similar to 

g f 


the one of interest, a reasonable guess of 


X 

s 


is readily available. 


o 

2. Housekeeping Power. This program option applies only to solar 
electric propulsion and is currently not integrated with the power degradation 
option. The program therefore does not allow the use of both options simulta- 
neously. The new program input -quantity controlling the housekeeping power 
simulation is DPOW, which is the ratio of housekeeping power to reference power 
and is given the symbol Ap: 


Ap = 


V p 


ref 


The reference power p is defined as the power input to the power conditioners* 
under reference conditions. (For solar electric propulsion, reference conditions 
are equivalent to having the total spacecraft in a perfect solar-simulator such 


*When the thrusters are at full -throttle. 
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that the simulated solar distance is one astronomical unit, the solar arrays are 
oriented normal to the source of energy, and no radiation damage or other de- 
gradation effects have occurred to any components). 

This spacecraft model assumes that the powers required by the operating 
components at each instant of time exactly match the power developed by the solar 
arrays: 


where p is the power developed by the arrays, p is the instantaneous power 

£L 

delivered to the power conditioners, and p^ is the housekeeping power, which 

is constant with time. Currently, all trajectories generated by the program 

using the housekeeping power simulation must satisfy the condition that p > p, , 

a h 

so that p remains positive; in other words, engine shutdown when p - 0 is 
not coded into the program. 

The power factor y is defined as the ratio of instantaneous power to 
reference power, 

n +4)/4 

y-p/p .= (1 + Ap) / a. d' -Ap 
ref 1 


where the summation quantity L is the power factor when p = 0. 

The spacecraft model has been expanded by deleting the old, total pro- 
pulsion system specific mass O' and replacing it with cn , the specific mass 

3 . 

of the solar arrays, and , the specific mass of the thruster subsystem. The 
quantity previously called the "propulsion system mass" is now expressed. 


m =(X p + a. 
ps a ao t 


ref 


<° t + °a & + AP,) P ref 


where p is p evaluated under reference conditions (i.e. , p =p + p, ), 
ao ao ref h 

and the total specific mass ot has been replaced by (a + a (1 + Ap)), also a 

t 8. 

constant since the housekeeping power option is intended to be used only in 
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simulations where the reference power is specified. The remainder of the 
analysis describing the optimization problem remains unchanged. 
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B. 


BOUNDARY AND TRANSVERSALITY CONDITIONS 


The basic boundary conditions and transversality conditions are described 
in Cl] ; the nomenclature used in the following analysis is also described in Cl] , 
except for the introduction of new symbols which are described in the text as they 
appear. 

In summary, when no intermediate-targets are present, individual trans- 
versality conditions are derived from the general equation 

° n 

d?r+ [/s • dV - A • dR + A dP + A ds+A df + A dg + A dc-h dt] =0 

Us T g c v o 

where it denotes the performance index which is assumed to be of the form 


IT - 17 (V , V , v , g, c, t , t ) 
“o 00 n n of 


such that 


d 77 — 57 dv +1T dv +17 dP + 17 dg + w dc + IT dt + IT dt 

v “ov °°n U n g c totn 
°°o “n n o n 

where ir denotes 5 77 /dx. The individual transversality conditions are obtained 
x 

by equating to zero the coefficients of all independent differentials. 

1. Deep Space Burn . In simulations of trajectories which are all-ballistic, 
the program is capable of simulating a single deep-space burn, or impulsive 
velocity-change, at any point prior to arrival at the primary target. The three 
components of the incremental velocity AV are independent variables of the 
boundary value problem, such that, at a specified time, the spacecraft velocity 
is incremented: 

*+ • - . 

R = R + AV 


The use of this program option is described in the Sample Problems and Results 
section under sample Case C, Multiple Ballistic Swingby Mission. 
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2. V„ Optimization, in LVI Mode, The optimization of the launch excess 
velocity' V a when using the Launch Vehicle Independent (LVI) mode of simulation 
is accomplished when the initial primer vector is forced, to vanish: 


A = 0. 

o 

This is accomplished by setting the values of these three independent variables of 
the boundary value problem to zero and turning their triggers off; the three com- 
ponents of the departure heliocentric velocity become independent variables in- 
stead, The program contains special logic to circumvent the numerical singularity 
associated with the null primer vector. 

3. Launch Asymptote Declination Optimization . This section contains the 
equations employed in the optimization of the launch asymptote declination. A more 
thorough exposition of this subject is found in [ 4 ] , in which some of the ramifica- 
tions of the theory are discussed. Historically, an incomplete form of this analysis 
was first published in the appendix of [l] . 

The assumed spacecraft and propulsion system models are as described in 
[l] and will not be repeated here. The launch vehicle payload capability is 
assumed to follow the simple exponential law 

m =b, e^cV -h , (20) 

O 1 o 


where b , b , and b are pre-determined constants for each launch vehicle and 
12 3 

v is a characteristic speed representative of the energy required to achieve a 
c 


specific escape trajectory. 


The velocity' penalty incurred with non-due-East launches from the ETR 
is shown graphically in Reference [5] as a function of the parking orbit inclina- 
tion. This velocity penalty Av. is adequately approximated with a quadratic 
curve fit of the form 

Av. = c i 2 + c„ i + c„ (21) 

1 1 2 3 
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Given a reference orbit inclination i and a circular orbit speed v , 

o 

the velocity penalty Av associated with a non-cop] anar departure from this 

g 

circular orbit to the desired hyperbolic excess velocity at a declination 6 is 

defined as follows,, Assuming the line of nodes of this reference orbit is an open 

variable, one may choose this variable to minimize the angle between the excess 

[6] 

velocity and the orbital plane. This minimum angle is 6 - i. Gunther L has 
shown that the minimum incremental velocity required to achieve a given v ro 
along an asymptote not lying in the orbital plane from a specified circular orbit 
is obtained from the solution to a quartic equation in the sine of the out-of -plane 


angle. Defining 


s = sin(6-i); P-v /v ; 

oo o 

2 2 

p = s (p +4) ; 

2 2 2 
q = s CL'S ) p ; 


x- [J(q/2) Z + (p/3) 3 + q/2_] 3 - [./(q/2) 2 + (p/S) 3 - q/ 2] 3 ; 

y = jp 2 / 4 - x ; (22) 

1 

w=| P /2 + y + J (p /2+y) 2 + 4^x/2 + Jx^7a + s 2 ^ ] , 

then Gunther’s solution for the magnitude of the minimum velocity impulse re- 
quired to accomplish the maneuver is 

v = Vq ./p 2 + 3-2 J (1 + pw-w 2 ) (2 + p w) , (23) 

and the penalty Av is the difference between v and the velocity increment re- 

g g 

quired if the out-of-plane angle were zero, i.e, , 


Av = v - ( „ 
a g V 


' v 2 + 2v ' 

00 r 
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Thus, the definition of the characteristic speed for those cases in which the 
asymptote declination lies outside the interval [— i, ij is 

v = /v 2 + 2v 2 + Av. + Av 
c V » o i g 

= v +v +Av.. (24) 

o g i 

The transversality conditions which optimize the launch phase are deve- 
loped as follows o The state and adjoint equations for the problem under considera- 
tion are precisely as formulated in [l] . The only difference in the optimality 
conditions is the format and content of certain of the transversality conditions . 

Specifically, these differences are due solely to the new definition of v which 

c 

is now a function of the direction of V as well as its magnitude „ Whereas be- 

fore the differential of v was, simply 

c 

dv = (v /v ) dv , 

C co c <» 

it follows from (24) that the equivalent formula now is 

dv = (gv /dv )dv + (dv /d6)d6 

C g oo oo g 

+ (dv /di+dAv /di)di 
g i 

where, from (21) 

dAv./di r: 2c i + c , 

1 1 u 

and, from (22) and (23), 

dv /d i = - dv /d 6 . 
g g 

The derivation of the partial derivatives dv /dv^ and dv /d6 is straight - 

g g 

forward although somewhat cumbersome. The equations for these partial 
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derivatives are as follows: 


2 2 2 2 

\v(3 + 2pv/-w )(c>p/9v c;) )+ (3p+2p w-3Pw -4 w)(^w/^v b ) 


3y v 

_£= -2_{p 

3 v v 3 v 




2 

+ pw-w )(2 + pw) 


}. 



2 2 2 
' ( 3p+2p w - 3pw - 4w) 

2v ./(! +pw-w^)(2 + pw) 


5w 
36 ’ 


where 


Sp/3v = l/v , 
^ ® o 



iw = l|i Z + f3x + ^ 3 y -] ; ; j 

36 2 l 3 6 U6 m: 2 \ 2 y / S 6 J 7 (Zw p/ ^ y) J ’ 

2,1 x /4 + s 


3 s/d 6 = cos (6- i), 


3 y/ 9v 00 = [(P/ 2 ) Sp/3v ro - Bx/dv^l^y , 


Sy/36 = - (3x/3 6)/2y, 


3_x _ 1_ 
&u 6 


(q/2)(8 9 /ou) . t( p/3i. g p /iuI + fiT / (q/2) 2 + (p/ 3) » +q/2 l 

/ftA? + G p/3) 3 3UjU J 


-2/3 


.-2/3 


_ I ja/g K^a/ &n ) + iP/- 3 L .C a .fi/ a u ) _ |SL ~([~y (q/2) 2 + (p/3) 3 - q/2 ] 
6L /^2/+(p/3) 3 
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with u = v or 6 , 

CO 

3q/o v co = 2 ps 2 ( 1 - s2 )(3p/ov co ) , 

3q/3 6 = 2p 2 s(l-2s 2 )(3s/3<5) , 

Bp/3 = 2ps 2 (op/By^), 

3p/35 - 2s(p 2 + 4)(3s/B6). 

This concludes the computation of these two partial derivatives. Some symbols 
used in this section of the document apply to this section onty; for example, p 
and q here are intermediate quantities in the analysis, and have no relation to 
power or degradation factor discussed elsewhere. 

The differential of g, the reference thrust acceleration, is 


dg = - (g/m Q ) dm Q 

where 

dm Q = - (b 1 /b 2 )e~ (V c //b 2 ) dv c , 
and the differential of V m may be written 


dV = 
00 


( v /v 

CO c 


) dv +(n xV )da+[(V 

v CO ' p co' 


x n )xV 

oo p' CO 


/|v 


x n 

CO p 


3 d6 


where n^ is a unit vector along the North Pole and a is the geocentric right 

ascension of V , and therefore the transversality conditions associated with 

V are as follows: 

00 


For optimum launch parking orbit inclination: 

f(3Av. /B i - 3v /3 6) = 0. 
i g 


(25) 
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For optimum launch excess speed: 


f(3v /Bv )-(A -V )/v = 0. (26) 

X gCO /V 0OO / CO 

For optimum launch asymptote declination: 

f(3 v /S6)-A • C(V xn )xV /|v xn | ] = 0. (27) 

g O'CCp'CO'OOp' ' ' 

For optimum launch asymptote right ascension: 

-A • (n x V ) = 0 . (28) 

o v p » 

In the above equations f is given by 


f=[k + k, - (1 +k) u.-gX /m ] dm /dv 
s t r f g o o c 

Due to the complexity of the equations defining °v /b 6 , the solution of equation 
(25) for the optimum value of i must be obtained using an iterative technique. 

The approach to the solution of the problem as formulated above differs 
in three basic respects from that of the problem where asymptote declination is 
ignored: (1) the condition (25) must be solved for the optimum parking orbit in- 
clination, given values of v and 6 ; (2) the asymptote declination 6 must be 

00 

introduced as an independent parameter and (27) added as an end condition of the 
problem; and (3) the evaluation of becomes somewhat more involved. The 
computation of V given v , A and 6, proceeds as follows: Denoteas e 

co co o 

the obliquity of the ecliptic such that the matrix 


$ = 0 cos e -sin c 


0 sin c cos C 
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operating on a vector expressed in ecliptic Cartesian coordinates yields the same 

vector in Earth equatorial coordinates,, Then the right ascension o; of the initial 

A 

primer A may be written 

a , = tan * [" (A cos e - A sin e)/A I , 

A L' yo zo xo J 


where A , A , A are the given ecliptic coordinates of A . Then, the right 
xo yo zo o 

ascension of the asymptote is set 


a= a, 


or a = + n 


and V is evaluated 

CD 



00 00 


cos a cos 6 
sin a cos 6 
sin 6 


4„ Extra-Ecliptic Missions, An additional set of boundary conditions has 
been added to the program for simulating extra-ecliptic miss ions „ 

The original set of boundary conditions is given on pp. 25-26 of reference 

/ 

[l] , which is now invoked by employing the program input quantity IOUT = 1 and 

corresponds to i, e, and r specified, Ci and o> optimized and f = 0, where 

P 

i is inclination to the ecliptic, e is eccentricity, r is perihelion distance, Q 

P 

is ascending node angle, co is argument of perihelion, and f is true anomaly, 
all evaluated at the final time. 


The new set of end conditions for extra-ecliptic missions is invoked by 
IOUT = 2 and is similar to the original set except that the final true anomaly is 
optimized rather than specified as zero, which corresponds to having the space- 
craft arrive at the optimum point of the extra-ecliptic orbit rather than being con- 
strained to arrive at perihelion. The specified final conditions are i, e, and a, 
where a is the final semi-major axis (in place of r ). Q and OJ are optimized 
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as before, using the transversality conditions 


s. 

D 

n 

zr i 

O 

n 

o 

(29) 

g 

n 

cn 

o 

O 
il . 

o 

(30) 


which are identical to equations (47) and (48) of reference [l] . C = R x A -R X A is 
the vector constant of the motion on a given trajectory segment, h is a unit vector 
along the angular momentum of the final orbit, and k is the z-axis unit vector. 
Optimizing f is accomplished by satisfying the- transversality condition asso- 
ciated with f: 


T (f) = ^ (A • R) + r 2 (A * R) = 0 (31) 

The program attempts to satisfy equations (29), (30), and (31) directly, since no 
closed-form solution has yet been found for casting these end conditions in the 

o 

form of AR and A.R similar to equations (54) and (55) of reference [l ] „ 

These end conditions possess poor convergence properties, a fact which the 
authors have experienced directly in attempting to obtain converged solutions 
using IOUT = 2. 


5. Multiple Target Missions - Errata . Two equations in the section, 
Extension of Formulation for Multiple-Target Missions of reference [l] are 
in error. Equation (116) of reference [l] should read 
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and equation (117) should read 


n-1 


977 


7 7 


m 3 m 
o o 

n-1 


= k +k +i - (1 + k ) v + (1 + k ) ; k 

s t J p 2'Q r n t -■ samp i 


i=l 


n-1 


r \ 

- k > k + g (1 + i k ) v -i k (1+ ' (k . - k , .)) 

tv- dropi x L' J t r n t t ^ samp i drop r ' 


i=l 


i=l 


c. 


AUXILIARY COMPUTATIONS 


This section presents equations employed in computations which are 
made after the iteration sequence involving the primary target is completed. 

1. Additi o nal Block Prin t Variables. A standard print block is employed 
for printing information at various points along a trajectory. Each standard 
block contains a total of forty parameters, which are described in [l] . 


The standard block may be augmented in two ways. When power degra- 
dation, as indicated by the input variable TPOWER, is simulated, a single line 
of information is automatically added to each block, as displayed in the Sample 
Problems and Results section, Case B, Comet Rendezvous Mission. When the 
input variable MPRINT is 2 or 3, three extra lines of information are generated 
per block, as displayed in the Sample Problems and Results section, Case C, 
Multiple Ballistic Swingby Mission. These two types of additionally printed lines 
may appear simultaneously. 

(a). Power Degradation. The single line of power degradation 
information contains eight parameters as follows: 


S 

LS 

DENSITY 

DPOWR 

DPOWD 

DEGRAD 

CHI 


Degradation time, s, since departure, in days, 
given by expression (4). 

Degradation time adjoint variable, A . 

s 

Density pai’ameter, d, used in expression (9) for 
the power ratio y, and given by expression (1), 
in AU -2 . 

3y/dr, in AU \ 

3y/dd. 

The degradation factor, q, given by expression (6). 

Solar array orientation angle X as given by, for 
example, expression (19), in degrees. 
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CHI REF 


Solar array orientation angle which the arrays 
would have if oriented for maximum power, in 
degrees. 


(b). Target-Relative Coordinates and Comet Magnitudes. The 
three extra lines which may appear via using MPRINT contain the following in- 
formation: 


R1 REL 
R2 REL 
R3 REL 

Cartesian components of current spacecraft 
position vector, with respect to the next astro- 
nomical body to be encountered along the 
trajectory in a moving coordinate system generated 
by that body, with the x-axis pointing outward along 
the body's heliocentric radius vector, the y-axis in 
the body's orbit plane in the direction of the body's 
motion, and the z-axis completing the right-handed 
orthogonal system, in kilometers, with the origin 
of coordinates at the body. 

VI REL 
V2 REL 
V3 REL 

Cartesian components of current spacecraft velocity 
vector, in kilometers/second, in the target -relative 
coordinate system described directly above (see 
R1 REL), 

RMAG REL 

Magnitude of R1 REL, R2 REL, R3 REL, in kilo- 
meters . 

VMAGREL 

Magnitude of VI REL, V2 REL, V3 REL, in kilo- 
meters/second. 

S/C NUC MAG 

Nuclear magnitude (of comet) of the next astro- 
nomical body to be encountered along the trajectory, 
as seen by the spacecraft. 


M xt = M + M„ log |r -R | + M log \R | 

N o 1 10 1 targ 1 2 b 10 targ 1 


+ . 03 cos' 1 [ Vs' (t Wl R> . 1 c o 

Kargl |B targ- R| 


where M , M 1 , and M 9 are magnitude constants 
(Continued on next page) 
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S/C NUC MAG 
(continued) 


S/C TOT MAG 


GEO NUC MAG . 
GEO TOT MAG 
ANG(V, R) 


ANG(V, XY) 


R1 REL ECL 
R2 REL ECL 
R3 REL ECL 

VI REL ECL 
V2 REL ECL 
V3 REL ECL 

RMAG ECL 


VMAG ECL 


associated with the target, and C° is the radians- 
to- degrees conversion factor e The arc-cosine term 
is the phase angle. 


Total magnitude of the next astronomical body to 
be encountered along the trajectory, as seen by the 
spacecraft. 


M - M 
T 3 


M 4 1OS 10 


|ft - R | + M, log, „ In | 
targ 5 10 targ 


where Mo, M^, and Mg are magnitude constants 
associated with the target. 


Same as S/C NUC MAG, except as seen by the Earth, 

Same as S/C TOT MAG, except as seen by. the Earth, 

Angle which (VI REL, V2 REL, V3 REL) makes with 
the positive x-axis in the target-relative coordinate 
system described under R1 REL, in degrees, 

ANG(V, R) = cos" 1 (VI REL/ VMAG REL) 

Angle which (VI REL, V2 REL, V3 REL) makes 
with the xy plane in the target -relative coordinate 
system described under R1 REL, in degrees. 

ANG(V, XY) = sin” 1 (V3 REL/ VMAG REL) 

Same as R1 REL, R2 REL, R3 REL except ex- 
pressed in the ecliptic coordinate system of 
date. 


Same as VI REL, V2 REL, V3 REL except ex- 
pressed in the ecliptic coordinate system of 
date. 

Magnitude of R1 REL ECL, R2 REL ECL, 
R3RELECL, in kilometers. 

Magnitude of VI REL ECL, V2 REL ECL, ' 

V3 REL ECL, in kilometers/second. 
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2. Swingby Continuation Analysis . Auxiliary computations are optionally 
provided, invoked by the NAMELIST input vector MOPT4, whereby ballistic 
swingbys past the primary target may be simulated. 

In one mode of program operation, invoked by MOPT4(l) > 0, single 
swingbys past the primary target may be simulated to up to ten post-swingby 
targets per case. 

In another mode of program operation, invoked by MOPT4(l) < 0, multiple 
swingbys along a single trajectory may be simulated, first swinging past the 
primary target and then subsequently swinging past more targets downstream 
along the trajectory. One multiple swingby trajectory may be simulated per 
case. 

In either mode of operation, the following basic assumptions are made. 

The swingby continuation computations are independent of the trajectory leg 
leading up to the swingby target, which may consist of an optimized electric 
propulsion trajectory segment (if the swingby target is the primary target), 
except that the arrival and arrival time at the swingby target are used in the 
determination of the swingby passage conditions. Each swingby maneuver is 
calculated under the assumption of the patched-conic approximation, and the 
swingby planet's sphere-of-influence is assumed to have zero radius as seen 
from interplanetary space and infinite radius as seen from the planetary vantage 
point. The passage time in the swingby planet's sphere-of-influence is neglected 
(taken to be zero in the heliocentric frame). 

Each swingby maneuver may be either unpowered or powered, and these 
two cases are discussed in the following sections. Since the unpowered swingby 
solutions are embedded in the wider class of powered-swingby solutions, tending 
to appear in pairs which are separated by a region of braking powered swingbys, 
the more general case of powered swingbys is discussed first. 
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(a). Powered Swingbys. This type of swingby maneuver is re- 
stricted to occur at the mutual perifoci of the approach and departure hyperbolic 
arcs; the powered phase is impulsive and the thrust is colinear (pro or con) to 
the velocity at closest approach. Whether the swingle is powered or unpowered, 
the trajectory segment leading up to the swingby planet has been pre-determ ined, 
this being the method by which the program has been designed to obtain swingby 
solutions. Therefore the swingby time and the arrival hyperbolic excess velocity 
V are known. In the following analysis, subscript A pertains to arrival at 
the swingby planet and subscript D pertains to departure. 

A basic assumption of the powered swingby problem posed here is that 
the flight time from the swingby planet to the next target is specified. This being 
so, the program is able to converge, by iteration, on some ballistic trajectory 
from the swingby planet to the next target having the specified transfer time, im- 
plying that the departure hyperbolic excess velocity at the swingby planet 

is thereby determined. Therefore, the heliocentric trajectory before and after ' 
the swingby planet is determined, and it then remains to perform the required 
computations pertaining to the hyperbolic arcs within the swingby planet's sphere 
of influence. 


and 


The closest approach distance is found by iteration as follows. Let 


oc = 1 + 

A 


r 

_E 


2 

v . 
coA 





+ 


r 

_£ 


v 


2 

k D 


A* 


where = |v | , v = [v | , r is the (unknown) passage distance, 

aJ p 

and /i is the swingby planet's gravitational parameter. Then the approach and 
departure hyperbolic bend angles are given by 
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—r— = cosec a A = sin (1 /o') 

2 A A 


D -1 . -1 

— ~ = cosec = sin (1/a; ) 


and these must sum up to the total bend angle, which is specified in terms of 

V . and V mT ^ : 

“A “D 


6 = -A + 

T 2 2 


- V_ A • V 


-1 °°A 

=. v 


v 

“A ro D 


Therefore, using r as the independent variable, the zero of the quantity 

P 

V • V 

-1 -1 . -i r cor) | 

F = sin (l/a a ) + sin (l/a D ) - cos j_ — - J 

00 A “D 


is obtained bj^ Newton’s iteration, using the derivative 


_9F = / -IN r «A^ g A + 

S r \ n t L /— 2 — : 

P VA _1 


2 / 
02 D 


X] 


Ad - 1 


When the iteration is converged, the passage distance r^ is in hand, and the 
impulsive velocity increment is computed, 


A v=/^ + v i 2 _im 
V r 00 D V r 
P 


where the square-root -quantities are the hyperbolic speeds at closest approach. 
The remaining parameters defining the planetocentric transfer are computed as 
follows. The inclination of the swingby orbit plane to the planet’s equator is given 


i = cos * (h • n ) 
' P 


where h is the unit vector along the angular momentum of the hyperbolic passage 

trajectory and n is a unit vector pointing toward the swingby planet’s north pole. 
P 
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The ascending node angle of the swingby orbit plane is computed as 

O = tan (-h /h ) 
x y 

and is placed in the proper quadrant by using the system library routine DATAN2. 
The argument of parifocus is given by 

— 1 * A 
co = cos ( i* ° r ) 
p n 

A A 

where r is the unit vector pointing toward the closest approach point and r^ 
is the unit vector lying along the line of nodes and pointing toward the ascending 
node. This is adjusted for the proper quadrant by the test, 

A A 

If h (r x r ) < 0, 0 ) ** 277 - 60 

z n p z 

In the right-handed planetary reference frame, the z-axis is toward the planet's 
north and the x-axis points toward the ascending node of the planet's equator on 
the ecliptic plane. 


(b)„ Unpowered Swingbys. This t 3 pe of swingby maneuver is 
considered to be a powered swingbj' having Av = 0. The program adjusts the 
post-SAvingby heliocentric trajectory segment, by iteration, until the swingby 
departure V magnitude equals the given arrival V magnitude. The primary 

03 CO 

independent variable in this iteration is the post-swingby transfer time to the 
specified target, which Avas held constant in the powered swingby case. Thus 
= v “ V oo * ant * the swingby passage distance is obtained from the 
formula, 



2 v 


CO 





The other orbital parameters are obtained from the same relations given above 
in the section, Powered Swingbys. 
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The program can generate multiple-revolution ballistic arcs, and a 
particular solution obtained by the program may not be unique, even for the 
same transfer time. All solutions are x’eachable, however, by means of in- 
putting. an appropriate initial velocity guess for the trajectory segment in 
question. 
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III. 


PROGRAM INPUT 


The following consists of a complete description of program inputs. 
With respect to the basic HILTOP report [l] , several input quantities have 
been added, some have been deleted, and some definitions have been modified. 
The input quantities which have been deleted are CALPHA, DIGIT, HALT, and 
WEIGHT, and so these names must no longer appear in NAMELIST MINPUT 
input data sets. CALPHA has been replaced by ALPHAA and ALPHAT. Also, 
Yi(4) and Yi(5) for i = 1, 2, 3, . . . , 70 have been deleted (See the description 
of "dependent parameters" in the section just ahead, Definitions of Input 
Parameters) and should therefore also be removed from input data sets. 
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A. NAMELIST 

Inputs to HILTOP are given through the NAMELIST feature of the IBM 
Fortran IV programming language. The input NAMELIST is named MINPUT, and 
every input required or used in the program is declared by name in the list. The 
general form for assigning an input value to a quantity is, simply, 

NAME = VALUE 

where NAME is the name assigned to the variable and is included in the NAME- 
LIST, and VALUE is a numerical or logical quantity consistent in form (i.e. , 
logical, integer, or real) with NAME. Unless otherwise specified, all MINPUT 
names commencing with one of the letters I through N represent integers, 
whereas all names commencing with one of the letters A through H or O 
through Z are double precision floating point numbers. Each NAMELIST case 
must begin with the characters 

& MINPUT 

commencing in card column 2 and followed by at least one blank, and end with 
the characters 

&END 

preceded by at least one blank. Card column 1 is ignored on all NAMELIST 
input cards. Multiple data assignments on a single card are permissible if 
separated by commas. Blanks in the variable field, VALUE, are taken as 
zeroes. A comma following the last VALUE on a card is optional on the IBM 
system. The order of the input data assignments is arbitrary; i.e., they need 
not be in the same order as listed in the NAMELIST. In fact, there is no re- 
quirement that any specific input parameter be represented in the input data 

set. If no value is included in the inputs for a particular parameter, the de- 

/ 

fault value is used (see Default Values). For other details regarding the 
NAMELIST feature, the reader is referred to the IBM System 360/Fortran IV 
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Language manual. NAMELIST cases may be stacked back-to-back indefinitely. 
A single NAMELIST error may wipe-out the remaining NAMELIST inputs. 
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B. DEFINITIONS OF INPUT PARA METERS 


Specific examples of the program inputs are given in the Sample Problems 
and Results section. Default-values of inputs are given in the next section. 

The program inputs, in alphabetical order, are: 


AAI 

Desired final extra-ecliptic inclination. Related to 
AE, AR, and IOUT. [deg] 

AE 

Desired final extra-ecliptic eccentricity. Related to 
AAI, AR, and IOUT. 

ALPHAA 

Specific mass of solar arrays, ct . Lkg/kw] 

a. 

ALPHAT 

Specific mass of thruster subsystem, o. . Ckg/kw] 

ALTITU 

This input variable is associated with program logic 
which has not been kept up-to-date, specifically, 
logic pertaining to optimum departure of a NERVA- 
type rocket from Earth orbit. This variable should 
be ignored. 

AN 

Trajectory-integration exponent in expression (37) of 
reference [l]. 

AR 

Desired final extra-ecliptic perihelion distance. Re- 
lated to AAI, AE, and IOUT. [AU] 

ASOL 

Array of five elements consisting of the solar power 
law coefficients a^ in expression (17) of reference 
[l] . ASOL(l) > 0 tells the program to use the input 
coefficients rather than the internal coefficients. The 
coefficients are normalized internally, and the program 
executes the iterations to produce the required remark- 
able points of the power curve (which are printed). 

BI 

Efficiency coefficient b in expression (15) of reference 
[lb Related to DI and El. 

B1 ' 
B2 , 
B3 

Launch vehicle coefficients b-,, b 2 , and bg in 
expression (2) of reference Cl j . Used only if MBOOST 
is negative, [kg, m/sec, kg] 
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CNI 

CNEX 

CSTR 

CTANK 

CTRET 

DI 

DMR.ETR 

DPOW 

ECI 

ECIX 

El 


Inclination to ecliptic of primary -target orbit. Input 
only when MO.PT3 = 11. Related to ECI, OMI, SAI, 

SOI, TPL EMUODD, and RADODD. [deg] 

Array of five elements, the first three of which may 
be currently used. Inclinations to ecliptic of inter- 
mediate-target orbits. Input CNIX(i) only when 
MOPTX(i) - 11. Related to ECIX, OMK, SAIX, 

SOIX, TPIX, EMUODX, and RADODX. [deg] 

Structural factor, k , in expression (7) of reference 

Cl]. s 

Propellant tankage factor, k in expression (6) of 
reference Cl ] . 

Retro tankage factor, k , in expression (10) of reference 

a], rt 

Efficiency coefficient d in expression (15) of reference 
[l]. Related to BI and El. Dan/sec] 

Retro engine mass, m rs , in expression (10) of reference 
[l]. [kg] 

Ratio of housekeeping power p^ to reference power P re f. 

The power transmitted to the propulsion system is that 

generated by the arrays less housekeeping power which is 

constant along the trajectory. The power output of the 

arrays normal to the sun at 1 AU is p + p, . 

J ref h 

Eccentricity of primary-target orbit. Must be less than 
unit}^ Input only when MOPT3 = 11. Related to CNI, OMI, 
SAI, SOI, TPI, EMUODD, and RADODD. 

Array of five elements, the first three of which may be 
currently used. Eccentricities of intermediate-target 
orbits. Input ECIX(i) only when MOPTX(i) = 11. Related 
to CNEX, OMIX, SAEX, SOEX, TPIX, EMUODX, and 
RADODX. 

Efficiency coefficient e in expression (15) of reference 
[l] . Related to BI and DI. 
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EMU ODD 
EMUODX 
GAMMAX 

GAP 


HOUR 

EBAL 

ENT PR 

IOUT 


Gravitational constant of primary-target. Input only 
when MO PT 3 = 11. Related to EC I, CRT, OMI, SAI, 

SOI, TPI, andRADODD. [m 3 /sec Z ] 

Array of five elements pertaining to the gravitational 
constants of inter mediate -tar gets. These inputs must 
be ignored at present. 

Maximum permissible value of the power function 7 
when MODE = 5. At solar distances less than the value 
for which 7 = GAMMAX, the solar arrays are assumed 
to be tilted such that 7 is maintained at the limiting 
value. . 

Propulsion-conier proximity tolerance-interval, Acr. 

See the discussion in the section Avoiding Corners in 
the Propulsion-time Function in reference Cl], When- 
ever the thrust switch function cr grazes the zero-axis 
within the tolerance |Acr | on any trajectory, an internal 
counter is incremented, and the trajectoi’ 3 ' is considered 
to be in the neighborhood of a propulsion-time corner. 
Positive value of GAP causes forced -thrusting case to 
be inserted, negative value causes bypass to next case, 
whenever the internal counter reaches the related input 
variable NHUNG. 

Hour-of-day of reference date (e.g. , 17.342D9). 

Related to MYEAR, MONTH, and MDAY. 

Ballistic option indicator. Setting IBAL 0 invokes 
option 1 discussed in the section Ballistic Trajectory 
Option of reference Cl] . 

Indicator which specifies print-length when the iteration 
in subroutine INTERP fails. Value of 0 causes short- 
print and 1 causes detailed-print. 

Extra-ecliptic mission indicator. IOUT = 1 or 2 
indicates that extra-ecliptic target conditions are 
desired, in which the iterator dependent variable 
triggers Yl(2) through Y6(2) are set equal to 1, 
and for which the input LAUNCH (which see) 
should probably be set to 1, and parameters related 
to LAUNCH also set appropriately. Ordinarily MOPT2 = 3. 
No retro stage may be employed. (Continued on next page). 
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IOUT 

(cont) 


IRK 


IRL 


IROT 


ISPIN 

ITF 


IT PR NT 


= 1 i, e, r specified; f = 0. 

= 2 i, e, a specified; f optimized. 

In the above, i = final extra-ecliptic inclination, 
e = final eccentricity, r = final perihelion distance, 
a = final semi-rnajor axis, and f = true anomaly at the 
final time. Final Q. and ce are optimized in both 
cases. Related to A E, AR, and A A I. 

Numerical integration option (currently not used). 


Primer-origin -proxi nitty step-size-control indicator. 
Value of zero causes bypass of control, leaving the 
step-size Au constant. See discussion in the section, 
Integration (Thrust) of reference [l ] . 

A non-zero value of IROT causes the input ecliptic 
projection of the primer vector and its time derivative 
to be rotated about the z-axis through an angle equal 
to the difference in longitudes of the spacecraft between 
the last trajectory of the previous case (or zero if no 
previous case) and the first trajectory of the current 
case. This feature permits one to use the initial ad- 
joint variables from a 2-dimensional trajectory as the 
initial -guess inputs for a 3-dimensional trajectory using 
the ephemeris option. IROT must be set back to zero to 
avoid undesirable rotations on subsequent cases. 

Spinner indicator. Not used at present. 

Provides normal termination conditions for rims which 
require more machine time than is estimated. The 
value specifies the number of machine-time seconds 
(CPU and I/O) required to execute the summary tra- 
jectory after halting the iteration-sequence, [sec] 

Does not apply if subroutine REMTIM is dummied. 

Indicator for special print from MINMX3 iterator. 
Non-zero value invokes print. 


JPP 


Jettison indicator j for electric propulsion system 

PS 

prior to primal-target retro-maneuver, as used in 
expression (8) of reference [l]. (Continued on next 
page). 
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JPP 

(cont) 


JPRINT 


JT 


KPART 


LAUNCH 


LOADX 


= 0 Propulsion system not jettisoned 
= 1 Propulsion system jettisoned prior to retro 
maneuver. 


Unit 11 printout-length indicator. A value of zero 
causes the iterator independent and dependent variables 
to be output only for each summary-trajectory; a value 
of one causes the same output additionall} 7 at each 
iteration of an iteration- sequence. 


Jettison indicator j for electric propulsion tankage 

prior to primary-target retro-maneuver, as used in 
expression (8) of reference [l]. 

= 0 Tankage not jettisoned 

= 1 Tankage jettisoned prior to retro-maneuver. 


Option for automatically selecting improved independent 
parameter perturbations for generating the iterator's 
partial derivative matrix. The option is invoked by 
setting KPART = N (N > 0), where N is the maximum 
number of allowed steps, as discussed in the section, 
Perturbation Step Size Selector of reference [l] . 
KPART must be set back to zero if not desired on 
subsequent cases. 


Launch mode selector, pertaining to the optimization 
of the departure asymptote declination, invoked by 
LAUNCH = 1. Related to X10, Y10, X17, and 
Y17. 

Intermediate-target initial-guess feature. Should be 
used with NSET(5) = 1, and then set to zero on the 
subsequent case. A non-zero value of LOADX will 
invoke this feature, whereby the primer A and its 
derivative A will be loaded into the iterator inde- 
pendent-variable arrays at each intermediate-target 
(Continued on next page). 
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LOADX 

(cont) 

MAXHAM 

MBOOST 


MDAY 

MODE 

MONTH 


provided that the trigger of the independent variable 
is on 0 The sole purpose of this capability is merely 
to generate an initial -guess for a multiple-target 
mission, where the values loaded into, the iterator 
arrays represent continuous A and A at each target. 

Maximum number of times that the program will print 
the warning message BAD HAMILTONIAN on any given 
computer run. 

Launch vehicle selector. 

= 0 ATLAS (SLV3X) /CENTAUR 

1 TITAN III C 

2 TITAN in C (1207) 

3 TITAN III X/CENTAUR 

4 TITAN JII X (1207) 

5 TITAN III X (1207)/ C ENTA UR 

6 SATURN IB/LM 

7 SATURN IB/CENTAUR 

8 SATURN IC/SIVB/ CENTAUR 

9 TITAN III X (1205)/ CENTAUR 

10 TITAN III B (CORE) /CENTAUR 

11 TITAN III D (1205)/ CENTAUR 

12 DELTA 

13 TITAN III D 

14 TITAN III D (1205)/CENTAUR/TE3G4 (2250) 

15 TITAN III E/CENTAUR 

16 SHUTTLE/TRANSTAGE 

17 SHUTTLE/DELTA 

18 SHUTTLE/AGENA 

19 SHUTTLE/CENTAUR 

20 SHUTTLE/CENTAUR/BURNER II (2300) 

NEG Use input booster coefficients b , b , and b 

JL ci O 

Day-of-month of reference date (e.g. , 26). Related to 
MYEAR, MONTH and HOUR. 

Power variation option selector. The value of MODE is 
equal to the option-number of the power -curve, discussed 
in the section. Electric Propulsion System of reference 
Cl 3 (which see). Possibly related to ASOL and GAMMAX. 
MODE = 1 has been eliminated. 

Month-of-year of reference date (e. g. , 8). Related to 
MYEAR, MDAY, and HOUR. 
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MOPT Ballistic option indicator. Using MOPT invokes option 2, 

discussed in the section, Ballistic Trajectorj' Option, 
of reference Cl], as follows: 

= 0 No action (use input A . A , and v ). 

' o o c o 

= 1 Generate ballistic solution with flyby end 
conditions. 

= 2 Generate ballistic solution with orbiter end 
conditions. 

Related to REVS. 

MOPTX Array of five elements, the first three of which may be 

currently used. This array specifies the target-number, 
or planet-number, of the successive intermediate-targets, 
and a value of zero indicates absence of the intermediate- 
target. A zero-entry must not precede a non-zero entry. 
Planet selection is the same as for MOPT2. MOPTX(l) 
pertains to iterator parameters X41-X50 and Y41-Y50; 
MOPTX (2) pertains to X51-X60 and Y51-Y60; and 
MOPTX(3) pertains to X61-X70 and Y61-Y70. Times 
at the targets are X48, X58, and X68. Not to be used 
unless MOPT2 ^ 0. 

MOPT2 Launch planet number and ephemeris-option indicator. 

= 0 Analytical planetary ephemeris is not used. 

^ 0 Analytical planetary ephemeris is used and 
the specific launch planet is selected as 
follows: (Continued on next page). 


. =1 

Mercury 

= 14 

Icarus (1978)* 

2 

Venus 

15 

Eros 

3 

Earth 

16 

Geographos (1983)* 

4 

Mars 

17 

Encke (1977)* 

5 

Jupiter 

18 

Encke (1984)* 

6 

Saturn 

19 

Encke (1987)* 

7 

Uranus 

20 

Halley 

8 

Neptune 

21 

Betulia 

9 

Pluto 

22 

Toro (1983)* 

10 

Ceres 

23 

Pallas 

11 

Input Target** 

24 

Juno 

12 

D’Arrest (1982)* 

25 

Vesta 

13 

Encke (1980)* 

26 

Astraea 


♦Year-value indicates apparition for which internal orbital elements are 
most accurate. 

♦♦Input corresponding orbital elements (see CNI, CNIX). None are available 
for the launch planet. 
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MOPT2 

(eont) 


= 27 Hebe 

28 Iris 

29 Flora 

30 Achilles 

31 Amor 

32 Hidalgo 

33 Alin da 

34 Grigg-Skjellerup (ID??)* 

35 Kopff 

36 Grigg-Skjellerup (1982)* 

37 Ganymed 

38 Ivar 

39 Beira 

40 Kepler 

41 Giacobini-Z inner (19S5)* 

42 Borrelly (1987)* 

43 Temps! II (1988)* 

44 Tempel II (1983)* 

45 Tuttle -G iacobini-Kresalc 

46 Schaiunasse 

47 Honda-Mrkos- Pajdusakova 

48 Giacobini-Z inner (1979)* 

49 Icarus. (3.987)* 

50 Toro (1987)* 

51 Geographos (1987)* 

*Year-value indicates apparition for which internal orbital elements are 
most accurate. 


MOPT3 Planet number of primary-target. Planet selection is the 

same as for MOPT2. If ephemeris is not used, MOPT3 
is used only for retro-stage mass computations, 

MOPT4 Array of ten elements, specifying up to ten post-swingby 

targets. Planet selection is the same as for MOPT2, and 
a value of zero indicates the absence of a post-swingby 
target. A negative value in MOPT4(l) selects multiple 
ballistic swingbys, rather than a set of single swingbys 
in which case also set MAXHAM = 0. Negative values 
(in absolute value) produce planet selection the same as 
for MOPT2. When MOPT4(l) < 0, the remaining elements 
of MOPT4(i) may be positive or negative. See the section, 
Swingby Continuation Analysis for details and Sample Case 
C for an example-case. Should be used only for primary- 
target flyby missions. Related to T2, MS WING, NSWING 
and XSWING. 
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MPGW 


Flag used in conjunction with the solar array degradation 
option. Value of zero results in the optimum orientation 
of the arrays relative to the sun line. A non-zero value 
forces the arrays to an orientation yielding the maximum 
power achievable at that instant. Related to TPOWER. 

MPRINT Indicator for printing the summary -trajectory (final 

trajectory of a case) as a function of time or for invoking 
extra printout. 

-0 Small-size block print at thrust switch points 
only (SWITCH POINT SUMMARY page). 

= 1 Same as = 0 , except expands to become a 
standard print-block of parameters for each 
computed point along the trajectory, including 
the trajectory extension controlled by the in- 
put variable TGO. 

= 2 Same as = 0, except each block contains extra 
lines consisting of target-relative coordinates 
* and target magnitudes. 

= 3 Combination of =1 and =2. 


MPUNCH Punched-card and trajectory-tape generation control. 


= 0 No special output. 

1 Punch final values of independent parameters. 

2 In addition, punch selected mission analysis 
parameters used for graphic documentation 
or other purposes. 


< 0 and > -100 Punch trajectory output used with the ASTEA 
program. The absolute value of MPUNCH 
determines the frequency of trajectory points 
output. 


s -101 Trajectory tape output used with the ASTEA 
program. The absolute value less 100 de- 
termines the frequency of trajectory points 
output. Related to NTAPE. 


MREAD 


MS WING 


MTMASS 


Card input option (iterator independent variables) 

= 0 No special cards input. 

1 The independent variables generated by 
a previous run by the MPUNCH option are 
input following the NAMELIST case, as 
discussed in the section, Program Output, 
of reference [l .1 . 

Array of ten elements, used only when running multiple- 
target ballistic swingbys, such that MSWING(i) corres- 
ponds to MOPT4(i) and selects the type of swingby 
maneuver desired at the respective swingby target. 

Used only if MOPT4(X) < 0. The shooting method 
(MINMX3 iterator) is used, and values of -1, -2, or 
-3 correspond to a swingby passage distance Initial 
guess of r = 63 (i.e. , continuous heliocentric velocity). 
Each element MSWING(i) may have any of the following 
values: 

= - 1 Go* directly for unpowered swingby; if and 

only if it fails, go for powered swingby having 
flight time T2(i) = initial guess. 

= - 2 Go directly for powered swingby only , having 
T2(i) - flight time, of post-swingby leg. 

= - 3 Go directly for unpowered swingby; then, 

whether it succeeds or not, go for powered 
swingby having T2(i) = flight time. 

= - 4 Go directly for unpowered swingby, but 
using initial velocity guess loaded into 
XSWING(j,i), j = 1,2,3, similar to 
MSWING (i) = - 1. 

= - 5 Same as = - 2, except use initial guess as 
in = - 4. 

*"Go for" means "attempt to obtain (solution)". 

Related to MOPT4, T2, XSWING, and NSWING. 

Mission-type selector pertaining to the primary target. 

= 0 Fljby mission 

1 Orbiter (high-thrust r etr o-m aneuve r without 
velocity loss). 

(Continued on next page). 
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MTIvIASS 

(cont) 


MUPDAT 

MYEAR 

NDIST 

NHUNG 

NORMAL 


= 2 Orb Iter (high-thrust retro-maneuver with 
velocitj 7 loss). 

3 Specified arrival excess speed v 

B n 

.If v =0, rendezvous mission 
»n 

If v > 0, controlled flyby mission 
00 n 

No retro-maneuver in either case. 

4 Orbiter (Electric propulsion system performs 
spiral maneuver. Arrival excess speed v 
must be specified as zero). 

Other parameters which may be related to MTMASS are 
DMRETR, CTRET, RPER, RAP, THRET, SPIRET, JPP, 
and JT. 

Flag indicating whether iterator independent variables 
at end of one case are to be used as first guesses for 
next case. 

= 0 Use input initial guesses. 

1 Use independent variables of previous case 
as first guesses. 

Year of reference date (e.g. , 19S2). Related to MONTH, 
MDAY, and HOUR. 

Identification number of celestial body to be used as the 
reference for the communication distance and angle 
measurement printed in the Extremum Point Summary 
Table. Identification code is the same as for MOPT2. 

Maximum number of propulsion-corner-proximity 
occurrences allowed in a given iteration-sequence. 

Related to GAP. 

Automatic adjoint-variable scaling. 

= 0 No action. 

1 All A and A are scaled such that 
becomes unity. 
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NPERF 


NPRINT 


NSET 


Identification number of end condition that is to be used as 
the performance index when employing the direct parameter 
optimisation feature (Improve Mode). The identification 
code is the same as the i in the Yi end condition array. 


Print selection flag. Permits selection of amount of 
printout desired on each case. 

= 0 Print only the case summary. 

1 Print switching point summary of final 
trajectory. 

2 Print MINPUT and case setup. 

4 Print trajectory summary on each iteration. 

8 Print partial derivative matrix each iteration. 

Combinations of options obtained by summing options desired 
If NPRINT > 15, printout consistent with NPRINT = 0 is 
obtained. If the sign of NPRINT is reversed to negative, 
the iterator independent and dependent variables addi- 
tionally are printed for every trajectory which HILTOP 
generates (including neighboring trajectories). 


Iteration-sequence control array. 


NSET(l) Not used for input. 

NSET (2) Not used for input. 

NSET (3) Maximum number of iterations permitted 
in attempting to satisfy (point) constraints 
in satisfy mode. If zero, no upper limit 
imposed. 

NSET (4) Flag indicating whether (point) constraints 

are to be satisfied prior to entering improve 
mode . 


= 0 Satisfy constraints first. 

1 Proceed immediately to improve mode. 

NSET (5) Maximum number of iterations permitted 
after entering improve mode. Setting 
NSET (5) = 1 causes iterator to be bypassed 
and computes single trajectory to obtain 
printout. 
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NS WING 


NSWPAR 


NTAPE 

OMI 

OMIX 

POWFIX 


Swingby continuation analysis option indicator. When 
zero or positive, NSWING is the number of full revolutions 
which the post-swingby heliocentric trajectory must have, 
and subroutine A CONIC (not MINMX3) is used to obtain the 
solution. When negative, NSWING has the same definition 
as MSWING (which see), and the MINMX3 iterator is used; 
NSWING must be used when MOPT4(l) > 0, and may be 
used when MOPT4 (1) <0. If MOPT4(l) < 0 and MSWING (i) 
= 0, then MSWING (i) will be set to the value of NSWING. 
Related to MSWING, MOPT4 and T2. 

Iterator independent-variable perturbation-increment 
control. 

= 0 No action. 

1 Allows the iterator to vary a given in- 
dependent-variable perturbation Ax 
whenever a neighboring trajectory is 
detected which has a different number 
of thrust switch points than the associated 
nominal trajectory. Ax is varied until 
the same number of switch points is 
achieved. 

Specifies the unit-number for the trajectory tape. 

Pertains to when MPUNCH * - 101. 

Ascending node angle (with respect to vernal equinox) 
of primary-target orbit. Input only when MOPT3 = 11. 
Related to CNI, ECI, SAI, SOI, TPI, EMU ODD, and 
RADODD. [deg] 

Array of five elements, the first three of which may 
be currently used. Ascending node angles of inter- 
mediate -target orbits. Input OMIX(i) only when 
MOPTX(i) = 11. Related to CNK, ECIX, SAK, 

SOIX, TPIX, EMUODX, and RADODX. [deg] 

Launch-vehicle-independent (i.e. , no launch vehicle) 
trajectory option in which the value of POWFIX is the 
spacecraft's reference power, [lew] 
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PSIGN • 


RADODD 

RADODX 

RAP 

REVS 

RPER 

SAI 

SAIX 

SOI 

SOIX 


Flag defining the sense of the launch hyperbolic excess 
velocity relative to the initial primer vector. A value of 
+1. results in the assignment of the geocentric right 
ascension cf the excess velocity equal to that of the initial 
primer vector. A value of -1. causes the geocentric 
right ascension of the excess velocity to be 180 degrees 
from that of the initial primer. 

Radius of primary target. Input only when MOPT3 = 11. 
Related to CNI, ECI, OMI, SAI, SOI, TPI, and 
EMUODD. [meters] 

Array of five elements pertaining to the radii of inter- 
mediate targets. These inputs are not used at present. 

Apoapse distance of capture orbit about primary target, 
[planet radii] 

Number of complete revolutions of the ballistic trajectory 
generated when the associated input MOPT is used. Must 
be a positive whole number. 

Periapse' distance of capture orbit about primary target, 
[planet radii] 

Semi-major axis of primary-target orbit (must be 
positive). Input only when MOPT3 = 11. Related to 
CNI, ECI, OMI, SOI, TPI, EMUODD, and RADODD. 

[AU] 

Array of five elements, the first three of which may be 
currently used. Semi-major axes of intermediate- 
target orbits (must be positive). Input SAEK(i) only 
when MOPTX(i) = 11. Related to CNK, ECEX, OMIX, 
SOIX, TPK, EMUODX, and RADODX. [AU] 

Argument of perihelion of primary-target orbit. Input 
only when MOPT3 = 11. Related to CNI, ECI, OMI, SAI, 
TPI, EMUODD, and RADODD. [deg] 

Array of five elements, the first three of which may be 
currently used. Arguments of perihelion of intermediate- 
target orbits. Input SOIX(i) only when MOPTX(i) = 11. 
Related to CNIX, ECEX, OMIX, SAIX, TPK, EMUODX, 
and RADODX. [deg] 
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SPIRET 


STATE 

STEP! 

STEP2 

TCOAST 

TDV 


TGO 


Retro-stage specif ic impulse (pertaining to the retro- 
maneuver at the primary target), [sec] 

Arra,y of six elements containing the Cartesian position 
and velocity components of the primary target. Use 
only when MOPT2 = 0 and the trigger settings of Yi(2) 
through Y6(2) are 0 or 1. [AU, AU/tauJ 
(tau = 58. 132440991 days) 

Thrust-phase computation step size, Au. 

Coast-phase computation step size, A/S. 

Array of twenty elements, consisting of the durations 
of the -coast phases corresponding to the coast-phase 
start-times input in the associated array TOFF. 

[days 3 

Time of occurrence of an impulsive deep space burn, 
in days from the start of the trajectory, which may be 
used only if the entire trajectory is ballistic (i.e. , electric 
propulsion is not permitted with this option, nor is a third 
intermediate target). Iterator independent variables 
X64, X65, and XG6 must be turned on, as these are used 
as the Av vector components of the deep space burn in 
EMOS. Also, set MAXHAM = 0. The following special 
feature is available regarding a first intermediate-target. 

If 1.D5 < TDV < 2.D5, then the burn occurs (TDV - 1.D5) 
days after passage of that target; if TDV > 2.D5, the 
burn occurs (TDV - 2.D5) days before passage of that 
target, [days] 

Ballistic trajectorj^-extension print option. When zero, 
no action. When positive, TGO = the number of days 
that the trajectory is to extend ballistically beyond the 
primary-target when no swingby-continuation is re- 
quested, and ballistically beyond the (last) post-swingby 
target when swingby-continuation is requested (in addition 
to the post-swingby trajectory segment itself). An3' nega- 
tive value will invoke printout of only the post-swingby 
trajectory segment or segments when swingby-continuation 
is requested. Applies also to trajectories with multiple 
swingbys . [days ] 
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TKRET 

TOFF 


TPI 


TPIX 


T POWER 


TSCALE 


T2 


XAKG1 


XANG2 


XSWING 


Retro-stage thrust, used only when MTMASS = 2. [lbs] 

Array of twenty elements, consisting of the times, in 
days from the start of the trajectory, at which imposed 
coast phases are to begin. Times must be in ascending 
order. Related to T COAST, [days j 

Time from reference date (MY EAR . etc.) to perihelion 
passage, for the primary target. Input only when 
MOPT3 = 11. Related to CNI, ECI, OMI, SAI, SOI, 
EMUODD, andRADODD. [days] 

Array of five elements, the first three of which may be 
currently used. Times from reference date (MYEAR, 
etc.) to perihelion passages, for the intermediate targets. 
Input TPDv(i) only when MOPTX(i) - 11. Related to 
CNEX, ECIX, OMIX, SAK, SOIX, EMUODX, and 
RADODX. [days] . 

Solar-cell degradation characteristic-time; nuclear 
electric propulsion radioactive-decay characteristic- 
time. Related to MPOW. [days] 

Iterator dependent -variable tolerance-interval scaling 
factor; scales all tolerances multiplicatively by the 
amount TSCALE. 

Array of ten elements consisting of initial estimates of 
swingby- continuation trajectory-segment flight-times , 
i.e., T2(i) correspond to MOPT4(i). [da} r s] 

Latitude of the launch site. Used only if LAUNCH is 
non-zero. Related to XANG2. [deg] 

Maximum parking orbit inclination, permitted by range 
safety considerations. Used only if LAUNCH is non- 
zero. Related to XANG1. [deg] 

Array of velocity vectors consisting of initial velocity 
guesses of a given post-swingby trajectory segment. 

Used only when either NSWING or MSWING has a value 
of -4 or -5. See especially the description of 
MSWING = -4. Velocity consists of exactly the same 
values as found in the VI, V2, V3 locations of the 
(Continued on next page). 
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XSWING trajectory block print (first block). Related to 

(cont) MSWING, NSWING, M0PT4, and T2. [AU/tau] 

XO Array of seven elements, the first six of which contain 

the Cartesian position and velocitj' components of the 
launch planet. The seventh element is not used for 
input. Used only when MOPT2 = 0. [AU, AU/tau] 


The following describes the iterator independent and dependent variable 
arrays of the boundary value problem. Input pertaining to the individual inde- 
pendent parameters is contained in the arraj^s XI through X70. The independent - 
parameter arrays have five elements for each variable, as follows (where i = 

1, 2, 3, ..., 70): 

Xi(l) Input value of parameter. Must be input regardless of 

trigger setting. If trigger is on (i.e., Xi(2) = 1), input 
value is used as initial guess of independent parameter 
and is varied at each subsequent iteration. If trigger is 
off, the parameter is not used as an independent para- 
meter and is not changed. 

Xi(2) Trigger indicating whether parameter is to be an in- 

dependent parameter in boundary value problem. 

Xi(2) = 0 Not an independent parameter. 

(Trigger is "off"). 

1 Use as independent parameter. 

(Trigger is "on"). 

Xi(3) Maximum change to parameter permitted in a single 

iteration. Should be a positive quantity. Used only 
if trigger is on. Units are same as that of the parameter. 

Xi(4) Perturbation increment used to compute partial de- 

rivatives by finite differences. Used only if trigger 
is on. Units are same as that of the parameter. 

Xi(5) Weighting factor. Should be a positive quantity. A 

value of 1. is generally recommended. The larger 
the weighting factor, the more the parameter is in- 
hibited from varying. Used only if trigger is on. 
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The independent variables are as follows: 


XI 

V> ■ 


X2 

A o (2) 

> Initial primer vector. 

X3 

A o (3) _ 


X4 

V> ' 


X5 

A o (2) 

> Initial primer derivative. 

X6 

\< S > . 


X7 

X 

VO 

Initial mass-ratio adjoint -variable. 

X8 

X 

r 

P r opuls ion-t ime adj o int - var iable . 

X9 


Not used. 

XI 0 

6 

Geocentric declination of launch hyperbolic 
excess velocity, [deg] 


There is no 
Xll 
XI 2 

XI 3 


conversion from input to internal units for any of the adjoint variables 

2 

Reference thrust acceleration, g. Cm/sec ] 
Electric propulsion system jet exhaust speed, 

c. [m/sec] 

Launch hjTperbolic excess speed, v 
[m/sec] 


X14 


I-fyperbolic excess speed at primary target, v 
Cm/sec] U 


XI 5 


Initial time, t Q , measured from the reference 
date (MYEAR, etc.), [days] 


X16 Time at the primary target, t , measured 

from the reference date (MYEAR, etc.), 
[days] 
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X17 


XI 8 
XI 9 
X20 


Launch parking oi'bit inclination, i. Used 
only if LAUNCH 0. Optimized internally 
by the program if both X17 and Y17 triggers 
are off. [deg] 


x 

o 



z 

o 


Initial spacecraft heliocentric velocity. Not 
required unless one of the three triggers is on. 
[AU/tau] (tau = 58.132440991 days) 


X21 


Constant thrust cone-angle, 0. Non-zero value 
invokes the constant-0 constraint. 0 < 0 ^180 . 
Zero-value implies that 0 is optimized along the 
trajectory (variable 0). [deg] 


X22 through X29 are currently not used (although some locations following X21 
are reserved for additional constant thrust cone-angles). 

X30 X Degradation-time adjoint-variable, 

s 

X31 through X40 are currently not used. X41 through X50 pertain to the first 
intermediate target, X51 through X60 pertain to the second intermediate target, 
and. X61 through X70 pertain to the third intermediate target. The correspond- 
ing intermediate- target parameters are ignored if the intermediate target is 
absent. Subscripts 1, 2, and 3 pertain to the first, second, and third inter- 
mediate targets, respectively. 


X41 

A x a) 

X42 

A x (2) 

X43 

A x (3) 

X44 

o 

A 1 (1) 

X45 

o 

^(2) 

X46 

• 

A x <3) 


Primer vector (at start of trajectory segment) 


Primer derivative (at start of trajectory segment) 


52 


X47 


X48 


X49 


X50 


Encounter speed at first intermediate target, 

v _ . Cm/see] 

“1 

Time at the first intermediate target, t^, 
measured from the reference date (MYEAR, 
etc.). [daj^s] 

Sample-mass factor, k „ , for sample- 

samp 1 

retrieval at first intermediate target. 

Drop-mass factor, k, , for instrument- 

drop 1 

package dropoff at first intermediate target. 


The independent variables X51 through XGO and X61 through X70 are 
identical to X41 through X50 except that they pertain to the second and third 
intermediate targets, respectively. A third intermediate target may not be 
present when simulating ballistic missions having a deep space burn (See TDV), 
in which case X64, X65, and XG6 are used as follows: 


X64 

Ax 

X65 

Ay 

X66 

Az 


Deep-space velocity-increment, [AU/tau] 


Inputs pertaining to the individual dependent parameters are contained in 
the arrays Y1 through Y70. The dependent-parameter arrays have three ele- - 
ments for each variable, as follows (where i = 1, 2, 3, ..., 70): 

Yi(l) Desired value of the dependent parameter. 

Yi(2) Trigger. If off (i. e. , equal to zero), the para- 

meter is ignored and is not considered a dependent 
parameter. Then the other two inputs pertaining 
to the parameter need not be input. If trigger is 
on (i. e. , not equal to zero), the parameter is 
considered to be a dependent parameter or con- 
straint. Certain of the parameters may have 
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Yi(2) 

(COJlt) 


up to three non-zero trigger settings. These 
will be discussed individually below. 


Yi(3) Tolerance of desired value (full interval width). 

It should be noted that the transversality conditions, which comprise some of the 
parameters, are developed under the assumption that all constraints are of the 
point constraint type. Therefore, the satisfy-mode is sufficient in solving any 
optimization problems for which a complete set of transversality conditions is 
available. 

The dependent -parameter •arraj’-s are as given below. T(x) represents 
"the transversality condition associated with x" and the function T(x) will have 
different values depending upon the constraints imposed on the problem. See 
NOMENCLATURE of reference [l] for definition of symbols and subscripts. 



Trigger 1 



Trigger 2 

Trigger 3 

Y1 

Ax [AU] 
n 

! 

r 

a [AU] 

Solar distance* [AU ] 

T(a) 

Y2 

Ay [AU] 
n 

t 

i 

e 

T(8 t )* 


T (8 t , 

Y3 

Az CAU] 
n 

T 

f 

i [deg] 



T(t n> > 

Y4 

Ax [AU/tau] 
n 

r 

f 

T(«) 

T(x n ) ■ 


T«) 

Y5 

Ay [AU/tau] 
n 

? 

? 

T(0>) 

T(y n ) 

optimal 

flyby 

V 

»0 

Y6 

Az [AU/tau] 
n 

T 

T(f) 

T <V 


T(X) 


♦Applicable only for two-dimensional motion in the xy plane. Also requires that 
MOPT2 = 0. 

Under Trigger 1 above, the first set of conditions applies to ordinary targeting 
conditions for position and velocity, and also to extra-ecliptic conditions to be 
satisfied when IOUT = 1; the second set of conditions applies to extra-ecliptic 
missions when IOUT = 2. T(f2), T(w), and T(f) are symbols for the transversality 
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NERVA 



conditions yielding optimum final node angle, argument of perihelion, and true 


an om al y , respect iv el y . 


Trig ger 1 Trigger 2 


Trigger 3 


Y7 

V 

n 

A 

vn 

m [kg] 

net 

Y8 

T (T ) 

r [days ] 


Y9 

Currently not used. 


Y10 

T(6) 

6 [deg] 

Used only if LAUNCH j* 0. 

Yll 

T(g) 

g [m/sec 2 ] 

I 1 

1 

lZj 

CD 

U 

Q. 

Y12 

T(c) 

c [m/sec] 


Y13 

T(v ) 

v [m/sec] 

coo 


Y14 


v [m/sec] 
“>n 

extra-ecliptic inclination [deg] 

Y15 

T(t o ) 

t [days] 
0 


Y16 

T(t n ) 

t [days] 
n 

t - t [days ] * 
n o 


♦Time transversality with flight time fixed is assigned to Y15 under trigger 1„ 
Y17 T(i) i [deg] , Used only if LAUNCH ^ 0. 


where i = parking 
orbit inclination. 


Y18 

T(x ) 

x [AU/tau] 


o 

o 

Y19 

T <V 

y [AU/tau] 
0 

Y20 

T(z ) 

z [AU/tau] 


o 

0 

Y21 

T(0) 

0 [deg] for 0 = 


constant with time. 
Y22 through Y29 are currently not used. 
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Y30 


(Degradation time) 


Trigger 1 Trigg er 2 

T(s) s [days] 

Y31 through Y40 are currently not used. Y41 through Y50 pertain 
to the first intermediate target. 


Y41 

Y4-2 

Y43 

Y44 

Y45 

Y46 

Y47 

Y48 

Y49 

Y50 


Ax [AUj 
Ay [AU] 

Az [AU] 

Ax 1 [AU/tau] T(x x ) | 

Ay x [AU/tau] | T(y ± ) 

Az ± [AU/tau] | T(z* ) 

v [ m / sec] 


optimal 

flyby 



Y51 through Y60 and Y61 through Y70 are identical to Y41 through 
Y50 except that they pertain to the second and third intermediate targets, 
respectively. 
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c. 


DEFAULT VALUES 


OF INPUT PARAMETERS 


The following is a complete, alphabetical list cf the default values of 
program input quantities having non-zero default values, except for the iterator 
arrays. All other inputs are zeroed. The default values of the iterator arrays 
Xi(l), Xi(2), Yi(l), andYi(2), for i ~ i, 2, 3, . 70, are zero, and the default 

values of Xi(3) through Xi(5) and Yi(3) for the same range of i are listed in 
the listing of program inputs of Sample Case C. 


ALPHAA 

15. 

NSET(3) 

300 

ALP1IAT 

15. 

NSET(5) 

300 

AN 

1.5 

NSWPAR 

1 

AR 

1 . 

NTAPE 

17 

BI 

.76 

POWFIX 

- 1 . 

CTANK 

.03 

PSIGN 

1 . 

CTRET 

1/9 

RADODD 

1 . 

DI 

13. 

RAP 

38. 

GAMMAX 

1 . 

RPER 

2. 

GAP 

.0001 

SAI 

1 . 

HOUR 

12. 

SPIRET 

300. 

IRK 

1 

STATE (1) 

1 . 

IRL 

1 

STATE (5) 

1 . 

ITF 

3 

STEP1 

.03125 

MAXI-I AM 

5 

STEP2 

. 125 

MDAY 

1 

TDV 

- 1 . ' 

MODE 

4 

TGO 

- 1 . 

MONTH 

1 

THRET 

400. 

MOPT3 

10 

TOFF 

20*~1. 

MUPDAT 

1 

T POWER 

10. **30 

MYEAR 

1975 

TSCALE 

1 . 

NDIST 

3 

T2(i) 

50* i 

NHUNG 

25 

X0(1) 

1 . 

NPRINT 

7 

X0(5) 

1.0000015 
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IV. 


SAMPLE PROBLEMS AMO RESULTS 


Three sample cases have been selected for presentation in the following 
pages. Each case was selected to display the use of one or more new features 
that have been added to the program. The first case is an extra-ecliptic mission 
that exhibits the effects of high launch asymptote declinations on launch vehicle per- 
formance and the low thrust trajectory. This case also displays the housekeeping 
power option that was recently added to the program. The second case is a comet 
rendezvous mission which includes the effects of solar array degradation due to 
radiation effects. The third and final case displays the HILTOP's powerful capa- 
bility for ballistic mission design and optimization. The specific mission chosen 
is a cometary flyby past Giacobini-Z inner followed by a deep space burn 10 days 
after passage, a return to and swingby of Earth (unpowered), a second swingby of 
Earth (powered), and finally encountering the comet Borrelly nearly 1023 days 
from launch. The tremendous flexibility for creating imaginative, multi-target 
mission profiles is demonstrated in this example. 
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A. 


EXTRA ECLIPTIC MISSION 


The objective of this mission is to place maximum pajdoad in a circular 
orbit of 1.001 AU radius inclined 45 degrees to the ecliptic. The mission is 
years (912.5 days) in duration and departs from Earth parking orbit on 
April 21, 1979. The solution given in this example is contained in the class of 
six-burn solutions which, for the mission duration assumed, tends to restrain 
the trajectory from deviating far from the nominal 1 AU solar distance. The 
specific case chosen uses the Titan III E/Centaur launch vehicle. The reference 
power (the power delivered to the power conditioners at 1 AU) is 20.35 kw 
and the specific impulse of the thrusters is 3000 seconds. The launch excess 
speed is optimized. The extra-ecliptic end conditions are invoked by setting 
IOUT = 1 and defining the desired values for AR, AE and AAI, the radius, 
eccentricity and inclination, respective!}'. 

This case exhibits the use of several optional features of the program. 

A total of 0.65 kw of power developed by the solar arrays is reserved for 
housekeeping (non-propuls ive) uses. This option is triggered with the input 
DPOW which is the ratio of housekeeping power to reference power. The power 
delivered to the power conditioners at distances below 1 AU is not permitted 
to exceed that delivered at 1 AU. This constraint is invoked by setting MODE 
equal to 5 and GAMMAX (the maximum permissible value of the power factor 
y) equal to 1. The effects of launch asymptote declination are included in the 
launch vehicle performance model by setting LAUNCH equal to 1. The equatorial 
inclination of the launch parking orbit is limited to a maximum of 36 degrees 
through the input parameter XANG2. Since the geocentric declination of the 
launch asymptote for extra ecliptic missions is usually much greater than this 
inclination limit, the solution will include a non-coplanar injection maneuver from 
the launch parking orbit. The declination of the launch asymptote is optimized. 
Finally, the option of inputting the coefficients of the power profile is illustrated. 
The inputs for this case are listed on the next page. 
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It should be noted chat the choice of final orbit radius of 1,001 AU rather 
than 1.0 AU was made to alleviate numerical difficulties arising as a result of 
the corner in the power curve at 1 AU. Neighboring trajectories terminating' 
on opposite sides of the corner point tend to possess different partial derivatives 
(i.e. , they will behave differently when subjected to the same perturbation). Con- 
sequently, if the final desired distance is exactly the point of the discontinuity, 
one might expect convergence retardation when the end conditions are nearly 
satisfied. 


AM INPUT Xl(2)«l.DO,X2<2W.nO / X3(2)=l.no / X4(2)‘»l.nO,X5(2>-l.DO 

X0{2)«l.n0;X7«l.D0 / X10(2)“l.n0 / Xll(2)*l.D0 / X12»2.9419D5D4 

X13(2) K l.D0,X15 s 3.103?0Dl,XI6=9.43537D2 

Yl(2)=l.DO/Y2(2)“l,no / Y3(2)=l.nO / Y4(2)=l.DO / Y5(2)=l.DO / Y6(2)=l.DO 

Y10(2)*l.D0 / Yll»20.35D0 / 3.n0 / Y13(2)»l.D0 

LAUNCH = 1 , M BOOS T= 1 5 , M TM ASS = 3,M00 E = 5 

MO PT2 = 3,M0 PT3 = 0,MY EAR= 1 9 79, M0 NTH = 3 ,MD AY = 2 1 

B I®. 6 3D 0,0 I = 0 . DO, CTANK = . 0550 0, 0AMMAX=1. 00,0 P0W = 3. 19 41031941030-2 
10UT-1, AAI =45.D0,AR=1. 00100, XAMG1®28. 500, XANG2=36 .00 
AS0L = 1. 4 3S 200, 0.0 0,-. 2 2 3500,0.00,-. 214 70 0 
Xl=-1, 8492690168360-02, X2 =-3. 444 5458696 870- 01, X3=-5. 292S946S082BD 00 
X4® 7.2838600516620-01, X5= 6. 4 899532772630-01, X6 =- 3 . 2 7 764 1 804 9 75D-0 1 
X10=-3, 80915072955 60 01,X11= 2 . 7905 74 06 12 560-04, XI 3= 5 . 32 86 06 534 64 5D 03 
& END 


MINIMUM 06NSITV = 0.0 MAXIMUM DPOwD = 1.A3820000000000000 00 

PROGRAM INPJTS 
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TRAJECTORIES WITHOUT PARTIAL OER I VA T I VES AND 3 TRAJECTORIES WITH PARTIAL DERIVATIVES REQUIRED FOR THIS CASE 








CM 

CM 

o 


44 

CM 

CM 

A 

44 

CM 

CM 

CM 

A 

44 

CM 

CM 

CM 

< 

CM 

CM 

CM 

CM 

A 

CM 

ui 






o 

o 

o 


© 

O 

O 

«“• 

o 

O 

O 

O 

- 

o 

O 

O 

O 


O 

O 

O 

O 

*• 

O 

o 

< 



N 



o 

o 

o 


o 

1 

o 

I 

O 

t 

o 

O 

O 

o 

O 

o 

o 

o 

O 

o 

O 

a 

© 

1 

O 

1 

O 

© 

o 

a 

o 


O 



♦ 


© 


. o 

o 

44 

© 

CM 

h* 

o 

44 

A 

CM 

A 

A 

A 

A 

CM 

o 

A 

A 

© 

CM 


© 


z 

UJ 


A 

♦ 

CM 


A 

N 

CM 

A 

© 

o 

N 

CM 

© 

© 

© 

© 

r- 

A 

N 

© 

© 

P* 

A 

N 


< 


u. 

X 


© 

♦ 

A 


N 

© 

A 

CM 

44 

© 

© 

A 


44 

A 

A 

© 

CM 

© 


A 

© 

© 

© 





44 


♦ 

44 

N 


© 


44 

O 

A 


c 

44 

© 

A 

O 

< 

*4 

A 

N 

A 

«• 

*4 

N 

S 


-I N 


X 

H 


# 

44 

O 


© 

CM 

44 

A 

A 

CM 

CM 

44 

CM 

A 

A 

A 

44 

44 

© 

5 

A 

*4 

A 

© 


U A 


u 



A 

44 

A 


© 

A 

44 

© 

A 

A 

A 

44 

PM 

A 

PM 

A 

44 


>- 

o> 

A 

44 

© 

N 


> => 



a 


© 

© 

© 


A 

© 

© 

© 

A 

*4 

© 

© 

A 

A 

O 

44 

© 

A 

© 

© 

*4 

© 

A 

© 


< a 

X 

4M 

o 

o 



44 

o 

CM 

K 

* 

© 

O 

44 

K 

4f 

A 

O 

© 

A 

* 

N 

© 

A 

A 

* 

A 

© 


a t 

< 

K 

a 

4 

• 

• 

* 

• 

• 

• 

• 

• 

0 

• 

• 

4 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


K N 

x 

A 

a 

o 

♦ 

© 

♦ 

o 

A 

< 

© 

A 

1 

A 

44 

♦ 

© 

1 

A 

CM 

A 

© 

1 

44 

CM 

A 

© 

CM 

44 






© 

O 

© 


44 

** 

44 

o 

O 

© 


44 

o 

O 

O 

44 


o 

O 

o 

44 

44 

© 

O 







o 

o 

o 

o 

o 

O 

o 

o 

o 

© 

o 

o 

o 

O 

O 

o 

O 

o 

O 

O 

© 

© 

© 

o 

o 







© 

O 

o 

Q 

o 

© 

O 

O 

o 

© 

o 

Q 

© 

o 

o 

o 

© 

O 

o 

o 

o 

© 

a 

o 

o 


o 


N 


A 

O 

o 

A 

A 

© 

A 

44 

O 

o 

© 

A 

44 

o 

© 

44 

A 

CM 

o 

44 

A 

A 

CM 

o 





O 


A 

o 

o 

© 

A 

. A 

A 

44 

O 

o 

© 

A 

44 

o 

A 

© 

* 

© 

o 

A 

© 

4 

o 

o 

© 


N 


z 


CM 

o 

o 

44 

A 

N 

A 

o 

O 

A 

© 

A 

o 

© 

A 

CM 

«r 

A 

o 

© 

K 

<3 

A 

o 

N 


< 


u. 


A 

o 

o 

A 

44 

A 

CM 

A 

o 

A 

CM 

CM 

A 

o 

A 

A 

«■ 

A 

o 

44 

A 

* 

A 

o 

CM 


a 




© 

o 

o 

© 

© 


4* 

O 

o 

© 

O 

44 

O 

o 

44 

O 

44 

A 

o 

N 

44 

44 

0 

o 

A 




a 


4f 

o 

o 

o 

A 

< 

44 

A 

o 

© 

A 

44 

A 

o 

<t 

A 

© 

44 

o 

A 

© 

© 

44 

o 

CM 

IS 

A 


UJ 

o 

o 

o 

o 

A 

© 

© 

s 

A 

o 

N 

CM 

N 

A 

o 

CM 

© 

A 

44 

o 

A 

CM 

A 

44 

o 

A 

< 

A 


k 

< 

o 

o 

o 

© 

A 

A 

A 

CM 

o 

O 

A 

A 

PM 

o 

44 

CM 


© 

o 

O 

N 


© 

o 

A 

X 

< 

N 

O 

X 

• 

• 

• 

• 

• 

• 

• 

• 

• 

4 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

* 

• 

• 


a. 

X 

J 

a 

> 

•* 

44 

44 

© 

© 

A 

© 

•* 

•* 


A 

© 

** 


44 

© 

A 




© 

A 



© 






CM 

44 

44 


O 

CM 

CM 

O 


o 

CM 

CM 

O 

44 

o 

44 


o 


o 

4 , 

cv 

o 

44 

O 


o o o o o 


o o O O o 


o o o o o 


o o o o o 


o 

z 

< 

A X 

o *• 

a. a 

w 

IS < ►- 

K a © x j 

< > J X u. 


z 

UJ ~ 

u r - 

unna i 

uj a j j a 

to 

X 

< 

a 

0 

3 

X 

1 < 

•* H 

X UJ 

HINNUr 

(0 <r j j »- 


z 

o 






UJ 

N 

V 




© 

A 

a: 




D 

D 

< 




N 

a 

X 



V 


X 

X 

UJ 


u 

© 

N 

D 

O 


o 

z 


if) 

O CM 

A 

-J 

o 

• 


z > 

-J 

u 

•J 

>- 

N 





a 

Z 





o 

44 





*- 

o 






& 





ui 


X 




•> 

X 

o 





o 

44 




or 


*- 



UJ 

»- 


< 



o 


* 

z 




IL 

A 

44 



*- 

O 


-J' 


UJ 




o 


z 

V- 

H 


Z 4« 


o 

< 

a 


-4 > 

-J 

o 

-I 

< 


o 

6 

O 

Q 

Q 


O 

O 

a 

o 

O 

o 

44 

44 

>4 



A 

<f 

44 

CM 

© 

o 

A 

© 

A 

A 


< 

A 

A 

© 

© 

o 

© 

44 

© 

© 


A 

rt i 

A 

A 

© 

o 

A 

© 

© 

A 


CM 

c 

>4 

(M 

o 

o 

A 

A 

A 

o 


A 

<y 

-4 

<f 

O 

o 

44 


44 

CM 


N 

CM 

N 

44 

CM 

o 

o 

© 

44 

44 


CM 

■o 

© 

© 


m 

A 

CM 

<» 

© 


A 

m 

A 

A 

A 

• 

• 

4 

• 

4 


4 

• 

4 

4 

• 

44 

44 

1 

A 

1 

© 

44 

1 


CM 

1 

A 

© 

A 

44 

44 

4< 

4 , 

CM 



44 




© 

© 

o 

o 

O 


© 

o 

O 

o 

O 

o 

© 

o 

Q 

© 


Q 

o 

o 

© 

n 

A 

N 

A 

•4 

© 


A 

© 

CM 

CM 

A 

«f 

© 


© 

< 



44 

A 

© 

© 

A 

CM 

A 

A 

44 


A 

A 

© 

44 


A 

© 

© 

© 

A 


A 


44 

A 

© 

<r 

© 

A 

© 

© 

U. 

© 

A 

© 

O 

N 

A 

44 

CM 

A 

44 

U. 

A 

■0 

o 

44 

© 

A 

o 

A 

44 

© 

O 

© 

A 

CM 

44 

© 

a 

A 

A 

N 



O 

CM 

A 

A 

© 

• 

4 

4 

• 

• 

1 - 




• 

4 

A 

N 

© 

44 

44 

A 


A 

© 

44 

© 


1 



1 

3 



I 


1 






cr 






© 

44 

44 

44 


X 

44 

44 

CM 

44 

4< 

o 

o 

o 

r» 


*- 

o 

o 

O 

o 

O 

A 

n 

A 

© 


T 

© 

o 

A 

o 

o 

© 

A 

© 

© 


O 

44 


N 

P- 

© 

A 

A 

00 

© 



A 


A 

o 

A 

N 

© 

44 

44 


44 

o 

o 

A 

A 

A 

CM 

o 

© 

N 


k 

A 

CM 

N 

© 

© 


A 

A 

CM 


A 

N 

A 

N 

P* 

O 

A 


A 

A 



© 

A 

44 

O 

CM 

CM 

© 

A 

O 



© 

CM 

© 



CM 

© 

o 

© 

o 


44 

< 

N 

N 

44 

• 

4 

4 

• 

• 


• 

• 

• 

• 

• 

© 

<r 


A 

o 


-* 

O 

< 

1 

A 

44 

1 



o 


44 


•4 

4* 

© 



© 


© 


o 


o 

o 

O 


o 

© 


© 


o 


o 

© 

© 


o 

* 


CM 


CM 


N 

© 

© 


o 

N 


A 


N 


44 

© 

44 


44 

A 


A 


© 


CM 

« 

© 


© 

O 


CM 


<r 


CM 

4? 

© 


© 

A 


CM 


A 


CM 

CM 

A 


© 

44 


CM 


44 


O 

44 

© 


A 

44 


© 


© 


A 

© 

44 > 



CM 

o 

44 

O 

N 


O 

N 

44 

O 

© 

• 

4 

• 

* 

• 


• 

• 

4 

• 

• 

44 

© 

A 

1 

o 

A 


44 

1 

44 

t 

O 

© 

44 

*4 

44 


44 



44 

CM 

CM 


© 

© 

o 


o 


o 

o 

O 

O 

O 

o 

o 

G 


© 


o 

o 

o 

o 

D 

A 

© 

© 




A 

CM 

<r 

A 

A 

N 

N 

N 


CM 


(M 

A 


A 

CM 

N 

O 

<t 


< 


A 

A 

A 

N 

O 

N 

© 

A 


44 


A 

<f 

© 

A 

o 

N 

44 

N 


A 


44 

A 

A 

© 


© 

44 

CM 




o 

A 


44 

A 

A 

A 



© 



(M 

© 

o 

N 

© 

44 

N 

o 

O 


CM 

O 

© 

A 

< 

• 

• 

• 

• 

• 


• 

• 

4 

• 


A 

A 

A 

o 

A 


© 

A 


N 

< 


O Q 
A 05 
03 K 
N <r 
<T 

CD A 
O CD 

*4 vrt 

A (M 


O O O 
h N O 
© © © 
A A A 

r> (\j o 
© <T ffl 
S •< O 
CM © w 
W if) o 


w 

3 

a 

x 


cm <o in © n 
l 

hmONh 

o o o O o 

I 

Q O O © © 

A A O' «T CM 

«* A <t CM fM 

n a) tfl n n 

n o w <f ifl 

O n rn o 

If) ffl O © O 

© O' m <t <t 

O W O N N 

• • • • • 

<r O' — ^ n 

I I 


r> o o o o 

i t 

0 0 0 0 0 

«« «n <o <n in 

«n O m K A 

O N O O O 

CONNO-- 
N — CM N © 

ifl S h N ^ 

O' C h N 

h N D O •« 

— M3 «-* O' - 

I 

'44 4* O *4 

0 o o o 

1 I 

o a o 

N N O 

h N O 

CM m CM 

N h S 

CM *i O' 

O -* cm 

a a o 

© p- o o a 


A 

n 

A 

o 

A 


*« H *4 O N 

I 


CM © 
O O 


0 O © 

1 t I 

o o a 

ANA 
CM © N 
n o »4 
A A CM 

P«OH 

o o * 
<T N O' 
CM O O' 

• • • 
O' a a 
i i 


a o 
a a 
a o* 

N p. 

n <r 
A A 

•4 H 

0 CM 

n a 

• 4 

n a 

1 l 


o o 
© n 
a <r 

A K 


© o 

CM A 
A 


A 

r> 

a 

x 


T 

u 


o 

<0 I 

A I 
A 

a o a 
N in n 
© A A 

c o n 
n a a 


p- 44 tf © A 

I 

H H H N ft| 

o © o © o 

I I 

o o o o o 

O A O' © A 

© O' o a as 

o -4 a A a 

CM O O' N o 

A CM O O* A 

P* N A O © 

o m k w o 

*4 O' © © -4 


CM * 
I 


0 o O O 

1 I 

0 O O o 
n o m o 

ID N K l\l 
© CM p>- a 
© A A © 
CO A <} *4 
A © N -4 
< A CM A 
A A © 4 
• • • • 
© A m *4 

1 I 


I 

O 


CM «4 O 

o o o 
I 

o o 

I CM A 
CM A CO 
AAA 
— © N 
O' CM A 
A © N 
© N O' 


O 

O' 

CM 

A 

O' 


» O N 
O' CM *4 © O' 


0 © 

1 o 


44 44 CM 

0 © o 

1 f 

o o o 

n <t n 

AAA 
<t A O 
<r © a 
t «4 o A 
• O a C 
I W C ffl 
» <t O O 

» 4 • • 

) © CM *4 


o o o o o 


o o 

A A 

o o 

A + 
A CM 
O' * 
O' © 

A -4 
© A 
• • 
O' CM 

I 


o o o 

4? cm CM 
O N ffl 
AAA 
<r O A 
<r A a 
AAA 
A A *4 
AAA 
• • • 
* O' * 


Z 

o 


A 

D 

a 

x 


x 

u 


I 44 44 44 CM 

. o o o o 

I I 

I o o o o 

4f 44 O' CM 

© CM N A 

> © © <r a 

O' 44 <r cm 

cm <r © A 

N N A A 

A A A *4 

I 44 A © «P 

• • • • 

N A A -4 

I 


r> o o r» o 

1 l 

o o o o o 

44 K O' O A 

«* cn A A o 

— A O N- O' 

A 4* N N CM 

tf N ON C 

A 44 o «f 44 

P» A < © A 

A A K A A 

• • • • ■ 

44 © ^ A 44 

I I 

(M -4 O 44 

0 o © o 

1 I 

o o o o 

44 <T <J 44 

W © N A 

AO© © 

44 A A 44 

© A O' CM 

A CM © 4« 

O O © A 

A © 44 o 4* 

© CM 44 o © 

I 


© o © 

I I I 
O O O 
r> A CM 
n a 4r 
CM cm A 

won 

- N n 

A ^ 

© A A 
CM A N 


© A A CM © 


O o 
t 

o o 

4 1 © 

A N 
A P- 
© 

O A 
A A 







44 

A 


44 

44 

4* 

*4 

© 

44 






0 

1 

O 


© 

© 

o 

O 

| 

O 

© 






1 

a 

o 


© 

© 

o 

o 

O 

o 






& 

© 


CM 

CM 

A 

< 

N 

© 






n 

o 


o 

© 

44 

N 

A 

© 




X 


A 

CM 


A 

44 

© 

A 

A 

A 




N 


© 

CM 


O 

A 

o 

A 

© 

N 




or 


CM 

O 


« • 

A 

o 

CM 

N 

© 

*4 



< 


© 

A 


A 

A 

•4 

N 

* 

44 


UJ 


UJ 


CM 

A 


O 

A 

44 

P- 

N 



X 

•4 


o 

© 

< 

o 

N 

O 

* 

© 

A 

o 

10 

M 44 

44 © IfJ 


• 

» 

• 

* 

4 

• 

• 

• 

• 

4 

<A 

N C 

J J 11 


o 

A 

4* 

o 

N 

0 

A 

A 

A 

* 

© 

V> 





1 



1 





1 


CM *4 4« © 44 
O © O O © 
i 1 

© o © © © 

©OAN4* 
r* A 44 A 44 

CM 44 44 O 44 

o o © © o 

1 1 

© © © O © 

44 « CM * N 

A A N *4 A 

CM 44 44 O 44 

o o © o o 

1 1 

© © © © © 

© © N cr © 

A CM 4* 44 P- 

*NAA4- 
P- A N © © 

ffl CM <M N k 

a) •< o c *• 

A A N N CM 

O © © ffl A 

44 © A N O 

K CM CM O A 

< A A * N 

© A A * A 

O P* A A © 

CM 44 < © A 

© CM <r N 44 

C O CM O C 

44 N © * © 

< 44 44 < O' 

A CM <f A <f 

cr A N O CD 

• •444 

^ © a k 

i 

4 4 4 4 4 
CM A CM ffl © 
1 1 1 

• 4 4 4 4 

CM N «0 A A 
. II I I 

o^^^JSuSw 

■OTP POOR QtJ&LS* 11 



67 


SWITCH THRUST OFF 


N 






w 

CJ 

w 

© 

PJ 

w 

W 

W 

< 

w 

w 

w 

N 

© 

N 

N 

OJ 

W 


N 

w 

N 

w 

© 

w 

tu 






o 

o 

o 


o 

o 

a 

o 

Y 

o 

o 

o 

O 


o 

O 

o 

o 


O 

O 

O 

o 


o 

< 




H 


o 

A 

o 

o 

o 

o 

o 

o 

o 

o 

o 

1 

o 

t 

o 

1 

c 

o 

o 

o 

O 


O 

o 

1 

o 

o 

1 

o 

o 

a 


U 


U 


ffi 

D 

<? 

© 

© 

a> 

ro 

c 

N 

© 

N 

ro 

O' 

ff 

N 

O' 

© 

O' 


N 

© 

44 

>4 

n 

«■ 



u 


z 


© 

© 

4 f 

o 

© 

© 

© 



O' 

«• 

<0 

to 

»P 

W 

© 

© 

r» 


W 

© 

N 

n 

© 

© 



< 


(L 

X 

© 


N 

© 

o 

N 

© 

K 

© 

o 

44 

© 

N 

*4 

«4 

© 

o 

N 


44 

44 


N 

N 

n 







* 

© 

44 

< 

o 

® 

o 

44 

© 

o 

© 

o 

«4 

o 

N 

44 


44 


N 

© 

© 

44 

© 

© 


j 

K 


X 

N 

o 

w 

44 

<r 

© 

N 

w 

44 

o 

o 

© 

44 

44 

* 

© 

w 

© 

44 


© 

N 

IT 

44 

CJ 

w 


cj 

© 


u 


•4 

© 

*4 

o 

10 


O 

44 

M 

ro 

ro 

© 

«4 

© 

© 

© 

w 

44 


© 

© 

N 

44 

w 

<y 


> 

D 


N 

fl. 

© 

N 

© 

cv 

© 

N 

N 

© 

w 

o 

*4 

© 

o 

O' 

N 

O' 

N 

© 


N 

© 

in 

© 

to 

© 


•< 

a 

T 

•4 

O 

ro 

*-« 

< 

pj 

< 

N 

4 * 

* 


* 

44 _ 



O' 

N 

to 

o 

* 

o 

N 

CD 

in 

* 

ro 

o 


a 

X 

< 

* 

cr 

• 

• 

• 

• 

• 

• 

• 

4 

• 

m 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

4 

« 



N 

N 

X 

© 

CL 


© 

© 

PJ 

1 

ro 

«* 

■© 

o 

1 

ro 

© 

N 

o 

1 

« 

© 

N 

© 

© 

4 

N 


© 

1 

© 







*4 


o 

o 

o 

44 


o 

o 

o 

44 

44 

o 

o 

O 

o 


o 


44 

44 

44 

© 

o 

o 







o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

© 

© 

© 

o 

© 

o 

o 

o 

o 

o 





© 

o 

O 

o 

o 

a 

© 

o 

a 

o 

o 

o 

o 

o 

o 

o 

o 

a 

a 

o 

o 

O 


Ct 


o 

N 

10 

tn 

© 

o 

ro 


Ui 

O' 

© 

© 

If 

w 

n 

o 

N 

© 

w 

n 

© 

© 

w 


© 

o 


•• 

U 

44 

w 

TO 

o 

o 

o 

w 

ro 

o 

n 

N 

© 

w 

o 

© 

44 

© 

w 

© 

© 

© 

o 

w 

o 


N 

z 

< 

N 

44 

o 

© 

44 

K 

44 

o 


m 

w 

c 

o 

© 

N 

w 

© 

w 

© 

o 

© 

44 

o 


«r 

u_ 

in 

44 

«r 

o 

o 

a 

44 

<T 

Cl 

© 

© 

<M 

© 

o 

N 

<T 

CJ 

© 

© 

<3 

N 

© 

© 

o 


a 


© 

© 

© 

o 

© 

<r 

© 

© 

o 

to 


© 

44 

o 

<f 

© 

© 

44 

© 

N 

© 


© 

o 



a 

© 

© 

N 

o 

o 

© 

© 


o 

o 

© 

W 

© 


W 

© 


© 


W 

© 


© 

o 

© 

in 

ID O 

in 

r- 

ro 

o 

w 

© 

N 

ro 

o 

<r 

44 

N 

ro 

o 

O 

O 

N 

ro 

© 

<3 

© 

© 

w 

o 

< 

V) 

* < 

p- 

fw 

r> 

o 

o 

© 

(V 

ro 

o 

© 

© 

rv 

44 

o 



W 

p>— 

o 

O 

© 

w 

n 

o 

X 

< N 

O X 

• 

• 

• 

• 

• 

• 

• 

• 

* 

• 


• 

• 

• 

• 

• 

• 

• 


• 

• 

• 

• 

• 


a i -i a > 


O' r- pj —t -* 


© N PJ 44 44 


© o n 44 4* 


■-< © JO © o 


o» © * »- *4 


vi 

o 

a. 


u 

d 

z 

< 


a 

13 

■*• v- 


o: n <5 i j 

< > j I u, 


Ui 

o 

3 


-J ID — 

U Zl- 

z «"• <r o « 

** > _J U -j 


u 

►4 

a 


w — 

u r- 

v ro n a x 


q o 

O © 
N © 
fv. lO 

© a 
<D 4 
© n 

K *4 

«r ro 


o o 

o o o 

o a a 

-* to o 
to n o 

O <4 <l> 
© P>. <» 

cvj cn 
© © © 
© to to 
in O' o 


ftl •< «f O N 

I I 

44 «-< O © 4 * 

o o o o © 

I 

o o o a o 

w w< O 21 S) 

W w ID O lO 

<} ^ n n io 

cd r- cj co <p 

© © N «? © 

o © in 


o © © 

t 1 

o o a 
o in w 
© o © 
© © o 
in © 

© n © 
««• n pi 
wo© 
o * © < 


© 

© 

in 

* 


o -« _» 

O Q O 
I ( 

o a o 

n r> © 

O S O 
n •* o 
o k ni 
© PJ <f 


Z - 


w - 
o o 

I 

Q O 
O N 

ro w 
<f pi 
«-* n 
© O' 

o w 

© <f 

© m 


w -« « © © 


^ o rg h 
o o o o o 
i 

o © © o n 
•*<©©© 
w © n n n 
<r <f o n n 

© O' © © p- 

© in © o in 

-4 O s N 


o o 

I 

O 0 
© © 
© N 
N © 

o •* 


O *4 o 

o o o 

o o a 
© © o 
© © to 

K O «P 
K. f> O 

<* © n 

<p © 

n ro © 
© © «♦■ 


© PJ 10 © *• 
I 

•* -4 O w ~ 

o o o o o 

o o © o o 

© n cd n o 

»o © ro © pj 

in o o »n n 

ro in *4 <y ro 

in in pj © n 

n w © © 


o o 
t I 

O O © 

© *4 W 
© © © 
© o © 
© < © 
(ON© 
— n in 
woo 


o 

© 

© 

© 

© 

n 

© 


O * « O N 


o o o 
( i 
o o o o 

© © © w 

© © 10 © 

© © w o 

<t © © 

W <0 N © 

© © © ro 

© »o © © 

© W N O © 


W -4 O 

o o o 

I 

© Q o 
O ro w 
o © o 

1 ft IO o 
lA <t © 
*-• 1-0 
© w -• 
A. ro 
44 © CO 
• • • 
~ n 
I I 


-I O o w 

o o o o o 

o o o a o 

© © © 4 O 

ro © ® o o 

in — k o o 

ro © ** © *>* 


0 o 

1 I 
Q O 
© w 
© ro 

o -* 
© © 
© o 
© N 
ro © 
© — 

• • 
© ffl 


© N <* 
W N 


<D © 
© 


ro 

© 

ro 

© 

w 

o 

n 

© 

«r 

PJ 

PJ 

o 

PJ 

© 

r0 

PJ 

© 

o 

w 

w 

ro 

© 

n 

© 

© 

44 

© 

4 P 


© 

p- 

© 

© 

o 


© 


© 

© 

© 


© 

pj 

N 

<P 

© 

• 

■ 

• 

• 

• 

N 

« 

• 

• 

• 

• 

H 

» 

• 

• 

• 

• 

H 

• 

• 

• 

• 

• 

to 

IO 

44 

© 

N 

© 

n 

© 

44 

•4 

r> 

© 

n 

w 

•4 

44 

© 

© 

ro 

© 

44 

ro 

*4 



1 


1 

D 



1 


l 



1 




3 


1 









a 






a 






a 






44 

44 

44 

44 

44 

X 

44 

44 

44 

PJ 

44 

X 

44 

4 « 

o 

4 * 

44 

r 

44 

44 

o 

44 

•4 

© 

o 

o 

r» 

o 

»- 

o 

o 

o 

o 

o 

►- 

o 

o 

o 

o 

o 

p- 

o 

o 

o 

o 

o 

© 

r> 

Q 

n 

ro 

X 

n 

n 

o 

CJ 

rv 

X 

n 

rv 

Q 

o 

© 

X 

rv 

1 

o 

n 

o 

A 

ro 

© 

© 

p- 

•-4 

u 

ro 

ro 

o 

ro 

p- 

cr 

<3 

44 

o 

44 

<n 

U 

<p 


pj 

44 

lA 

© 

N 

PJ 

lO 

o 

»- 

© 

p» 

© 

44 

«p 

»- 

•/l 

© 

K. 

© 

© 

W 

© 

<r 

ro 

m 

PI 

© 


«n 

© 

<r 

44 

© 

o 

ro 

© 

44 - 

44 

o 

<r 

ro 

© 


44 

o 

N 

P. 

_ 

W 

© 


44 

44 


31 

© 

N 

44 

N 

PJ 

* 

© 

O 

PJ 

© 

N 

* 

© 

ro 

© 

w 

© 

© 

© 

N 

© 

© 

© 

© 

© 

© 

44 

© 

© 

© 

n 

O 


N 

© 

© 


ro 

© 

© 

w 

N 

ro 

N 

ro 


w 

PJ 

ro 

© 

© 


to 

© 

© 

44 

O 


IO 

© 

ro 

w 

O 

44 

© 

W 

W 



44 

N 

© 

o 

© 


© 

© 

o 

c 

n 


© 

ro 

44 

ro 

ro 

w 

•«r 

© 

© 

o 


PJ 

* 

N 

o 

o 




w 

© 

© 


© 

44 

44 

ro 

© 

• 

• 

• 

• 

• 


■ 

• 

• 

• 

• 


• 

• 

• 

• 

• 


• 

• 

• 

• 

• 

PJ 

© 

<* 

© 

PJ 

1 


w 

© 

© 

•• 

P) 

1 


w 

© 

1 

7 

N 

pj 


ro 

© 

1 

( 

p- 

PJ 

w 

44 

o 


44 


PJ 

44 

© 


44 


PJ 


o 


44 


pj 


© 


*4 


o o © 

! » 

o o o 
© © © 
© ro o 
©NO 
< o * 
WO© 
© © « 
O N 44 

© < N . 


O 

© 


W *4 © O -4 

o o o o o 

O o O O Q 

<r © w n n 

N © © © <t 

pj ♦ © © 

O <n O <u n 
n o> ® 

© *4 © K) x 

© © * o o 

ro © © o in 

• • * • » 

W w ro — © 

t \ 

•4 O W *4 

o o o 


© o 

I 

o o 

N © 
© N 
© *r 
© © 


o o o 

N w © 
•- N w 
W < *-. 

o © ro 
O P» O' 
N © ® 
© r- C\j 
O' «T N 


»4 #0 

I 


o o 
1 

O ro 

— N 
<♦ o* 

— o 
© O' 
N W 
O' ro 
© ~ 
© ro 


o w 

o o o 

oori 

£ft 'ft ftl 

N «-• r> 
© ru ro 
O' < 

O ro © 
own 
<V N *4 
W O *4 
• • • 
•• -4 n 
I 


w — 
o o 
« I 

o o 

N -< 

© 4 f 

© o 
® N 
W O' 
W <f 
N © 

ro ^ 


O 44 

o o 

o o 

© O' 

o* © 

w w 

© 0* 

© © 

w o* 

r0 O' 

© O C\l 


id a 

J JQ. 

^ IO 
1 

44 

1 

o 

N 

< 

10 

i 


© 

N 

PJ 



O 

© 

w 

4 

44 

1 

O 

© 

44 

© 

1 

PJ 

o 

© 

© 


44 44 

44 

44 

•4 

44 


44 



44 

44 



44 

44 


44 


44 

o 

44 

*4 

o 


X 


o o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

© 

o 

o 

© 



o o 

o 

a 

o 

D 

a 

o 

o 

o 

o 

o 

o 

o 

o 

a 

a 

o 

o 

a 

o 

o 

o 

n 

o 

a 


n id 


o 

o 

n 

© 

N 

© 

© 

ro 

PJ 

N 

c 

o 

ro 


o 



© 

N 

(O 

© 


o 


♦t o 

© 

44 

44 

< 

© 

IO 

44 

PJ 

© 

© 

N 

PJ 

© 

© 

© 

<r 

PJ 

n 

N 



<r 

© 



PJ 4 T 

—4 

PJ 

44 

w 

<3 


PJ 

N 

© 

n 

O 

N 


© 

o 

© 

N 


© 


© 

o 




© © 

o 

>o 

o 

© 

PJ 

44 

© 

44 

o 

© 

N 

© 


o 

o 

© 

© 

© 


© 


N 

© 

X 


N © 

© 

© 

o 

N 

n 

o 

© 

© 

in 

w 


«f 

© 

© 


© 

<r 

© 

o 

W 

o 

N 

<r 

» 

< 

© 44 


o 

ro 

© 

o 


© 

© 

o 

© 

© 

N 

n 



N 

N 

© 

© 

O 



© 

** 

H 

W O 

r> 

© 

o 

PJ 

o 

© 

© 

<3 

© 

r> 

PJ 

C 

© 

© 

® 

© 


© 

o 


© 

O 

w 

X 

ID 

© © 

44 

ro 

© 

© 

© 

© 

n 

ro 

© 

© 

44 

44 

© 

© 

© 

© 

44 

© 



© 

O 

© 

UJ PJ 

W U X 

• • 

• 

• 

• 

• 

• 

• 

% 

• 

• 

4 

• 

« 

* 

• 

• 

• 

4 

• 

4 


4 



© or 

J J H 

© o 

w 

* 

l 

«r 

© 

* 

1 

N 

1 

* 

i 

W 

1 

© 

® 

© 

1 

N 

1 

© 

© 

ro 

© 

N 

1 


44 

« 

1 

© 

1 

n 



W *4 

44 

44 


PJ 


44 


44 

w 

44 



44 

w 

44 

44 

44 


PJ 


-4 

44 

44 



© © 

o 

© 

o 

O 

© 

o 

o 

o 

o 

o 

o 

O 

a 

o 

o 

o 

o 

o 

o 

© 

© 

o 

o 



o o 

D 

a 

o 

o 

a 

a 

o 

o 

o 

o 

o 

o 

o 

a 

o 

o 

o 

o 

o 

o 

o 

o 

o 



n © 

© 




m 


© 


PJ 


© 

PJ 

© 

« 


PJ 



© 


PJ 





© ro 

© 

ro 

ro 

PJ 

© 


n 

® 

© 

© 

© 

© 

© 

ro 

o 

© 

© 

o 

© 

N 

© 

c 

N 



o •< 

© 

pj 

© 

PJ 

© 

© 

pj 

N 

PJ 

N 

N 

N 

ro 

N 

PJ 

* 

N 

W 

© 

© 

© 

ro 

© 



N <3 

N 

w 


© 


pj 

w 


PJ 

© 

4 * 

O 

o 

N 

© 

44 

o 

& 

o 

4 t 

© 

© 

© 



O N 

PJ 

4-4 

ro 

© 

n 

© 

4 . 

44 

© 

© 

© 

© 

N 

S 


U) 

© 


CV 


© 

© 

N 



O O 


© 

© 

© 

<3 

4 f 

© 

PJ 

44 

© 

N 

© 

44 

N 


© 

© 


© 

IO 

•4 

© 

O 

ID 


44 © 

© 

© 

o 

© 

PJ 


© 

© 

© 

© 


< 

PJ 

N 

w 

© 



© 

N 

© 

n 

n 

X 


o « 

ro 

© 

© 

n 

N 

o 

© 

© 

© 

in 

N 

4 < 

o 

44 


© 



* 

O 


© 

.4 

P 4 4 f 

o © 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


• 

• 

4 


• 

• 

• 

• 

• 

I- at 

J J (L 

4 

2 

PJ 

44 

« 


N 

<0 

•4 

N 

in 

44 

PJ 

PJ 

© 

© 

« 

N 

PJ 

© 

N 

4 ? 

© 

w 

© 





1 

i 




1 



\ 


1 



1 

1 

1 

1 



1 

l 

1 


68 


n 


N N N O « CM *< «-* <V 

m O O O -4 O OOOOO 

© 1 I I <11 

< >- ft ft O ft O c a o o o 

g, u y © n cn <? <? 

y X ui N ft O © O <f © © © 

< IL* O « N lO © © N © ® 

Ml © © — «f O H ffl H O' N 

J K X ►- © «* — tv O O - « Ifl 

IJ {/) u -« O & N D « O 

>3 f-a # in m\i c o — c tr> «r 

<K x mO ff in « n o < © * n o< 

c: z < * tf ••••• ••••• 

^ K Z </> 0- ® ® © CM © © •- © © N 


O f* 

. - u 

•- T 

< u. 

ec 

a 

O CO Ui © 

« tf) * C 

x < r- o x 

c x y a > 


o 4 o 4 O O O 

o o o o o 
< I 

o a o o o 
© ^ in o in 
NOWON 
# m « o 0 
n in o o in 
•- <y O' o o 

© <f © O © 

<t © cm o o 
o w cl o o 

• • • • • 

© iT> <r o* o« 


o »-> o •-* «•* 

O o o o o 

< I I 
ft O ft O O 

ft n r> <■ \n 

>o in n in c 

<t c o a>- m 
N S O m O 

© <J o. C 1/3 

OONOn 
O W ID O O 

O « K O' O' 

• • • • • 

** + 1/3 © © 





w •* O o -4 


04 o4 O O ft 

ui 


o o a o o 


o o o o o 

1 

y 


O O O Q O 


ft o ft a o 

z 


o n n n w 


O © ft <t © 

< 


cd ft m 


k © <n © © 



m in «v co © 


K ft O CM CM 

(0 T 


-4 o X) CM O 


04 CD 04 ft ft 

O >- 


N (O N O' O' 


N. CD © © © 

a a 


o 4 <0 o n o 


© © © O <* 

y 


s rt d o o 


s ai o a O 

O < K 


O' c* w o m 


CM o4 O © 

tr ft © x j 


4 • 4 • • 


• 4 4 4 * 

< > J X U. 


N cm n 04 04 


< 1/3 CM 04 04 

.1 1 



04 o4 o (Vi 04 


04 ft O CM -4 



o O o O o 


O o O O O 



OOOQQ 

z 

ft O ft ft o 



K o <* O' o 4 

0 

© \f) * © © 

ttl 


© ft o h- © 


© © ft ^ © 

o 


f> O' N (V <f 

►- 

N © ft © © 

D 


© r> n o 

w> 

ft O O © © 

H 


o o w ® <j 

D 

CM N © ft © 

if ^ 

2 

© © tn © * 

a 

© © O 04 © 

tu U 13 

o 

i/l N o in m 

r 

K © CM CM O 

ft o z 


® « ■< o m 

»- 

© *4 © O © 

ONV1 JO 

K 

4 0 4 4 4 


4 4 4 4 4 

z > J VJ J 

tn 

ft © CM 04 04 
1 1 

• 

n © 04 cm © 


a 


s 



X 

04 04 O C\J 04 

o 

04 04 © 04 04 


K- 

O O O ft o 

to 

O O o O r> 



1 

U 

1 

2 ' 

X 

ft © © © © 

Ui 

ft r, a ft o 

o 

u 

•4 O N N n 

-1 

© O' .4 CM © 

M 

' H* 

4» o 4f rt © 

< 

n K f> k in 

►» lli 

04 

H N <J <f H 

ff 

<D irt N 

< o 

* 

© K © 04 © 

♦- 

N N o CM 

Z 3 

to 

f«- ft CM © CM 


© © ft © © 

*- »- 


© © cm <r> m 

u. 

© ft ft o o 

J LU ~ 


o r> O m o 4 

o 

© <t © CM N 

u z »- 


© h- © 04 cm 


^ © CM ft © 

Z -4 « O < 


4 4*44 

ft 

• • 4 • 4 

* i J u J 


cn 04 o 4 n 

1 

z 

tu 

<MO - N W 
1 1 



CM 04 O 04 


© 04 O 04 



O O © O 


© o o o 

> 


1 1 


1 1 

h 


O ft o o 


ft O O ft 



mv n •* 


© o o o 

u 


to < <f * 


ft © o © 

04 


© © © © 


© N ft © 

a 


© 04 © 04 


CM N © o4 

»- 


CM CO ft CM 


© © © © 

z 


W o4 N © 


cm r* © ♦ 

Ui *4 


W O CM 04 


© O CM CO 

U T - 


ft CM © O © 


4 t © N O N 

u ft ft a x 


4 4 * 4 * 


4 4 4 4 4 

me J jo. 


H #) 04 O CO 
1 


© * CM O © 

to 

04 


O O O 04 


O 04 o O 04 

K 


© o o o o 


o o o © o 

< 


ft O ft ft ft 


1 

ft ft ft ft o 

<X 


© 04 CM © CM 


© * O N CM 

O 


s .4 <o < n 


ft © 04 © © . 

•> 


©©•40 © 


© <f a> © © 

< 


♦ tn ft n <T 


© © © © h* 

X 


© <M ^ N ft 


O CM <■ ft CM 

1 < 


© 04 © m ft 


04 © <y <\J CM 

04 *- 


O CM CM © © 


O <t K — © 

x tu 


o s < o © 


O CM © l/> © 

UKVNUX 


• • • • • 


• 4 4 • • 

to a j j »- 


04 CM 04 ~ © 

i I i T 


— © 04 04 © 

1 1 



(M 04 o4 04 o4 


CM 04 o 04 04 



o o o o o 
1 1 

H 

O O © O O 
1 



ft ft ft O ft 

UJ 

ft O O ft ft 



© © f- 4f CM 

o 

o © © — N 



© 4 * tn © 

a 

O O <t © < 



4* © K ft © 

< 

O ft o. O o. 



N < © 


O — -4 K <J 



©to © © © 


o © <f cm 



© ft to to © 

¥- 

© N © O -• 

tu 


© N ft ft ^ 

0 

NO-MI 

z « 


o © © © © 

a 

o4 o CM « ft 

04 *4 04 O <0 


• •444 

z 

• • • 4 • 

a j j a. 


© K N W N 

04 

© ft 04 ft © 


I I 



69 


CASE 1 iterator summary 



m4 


** 


** 

** 



*■* 


** 


*■* 

o 

o 

© 












o 

O 

o 












© 

* 

o 












n 

ft 

o 












© 

© 

o 












© 

n 

ft 












(y 

© 

© 












tf) 

© 

o 












n 

ft 

pp 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

• 

• 

• 

• 

• 







• 



© 

ro 

ro 

o 

o 

© 

o 

o 

o 

o 

o 

o 

o 

o 

CM 

Z 


to 

© 

a 

in 

o 

< 

a 

c 

a 

o 

a 

ft 

J 

tu 

o 

ft 

o 

— 

mm 

1 

o 

1 

o 

1 

o 

o 

u 

i 

-» 

a 

u 

X 

M 

O 

cc 

(U 

A" 

Ci 

a 

o 

LU 

«-• 

UJ 

r 

o 

a 

X 

o 

a 

O 

o 

o 

o 

a 

O 

ft 

> 

ft 

■J 


a 

a 

v. 

a 

X 

a 

x 





• 










to 

O 

tn 

o 

in 

o 

in 

o 

in 

o 

in 

o 

in 

o 



*M 

fJ 

0 J 

ro 

ro 

<t 

* 

in 

in 

© 

© 

ft 


to 

ft 


-* 


- 


** 



•* 






o 













a 
















o 

o 













o 

© 

© 














ft 

© 














ft 

© 













ft 


ft 














ft ft 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 


» 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

N 

to 

o 

o 

o 

o 

o 

o 

o 

o 

© 

o 

o 

o 


w 



w 


w 


w 

W 

w 

w 



w 

* 

© 

z 














> 

mi 















u 














ft 

tu 














mi 

Q 














tu 

• 













7 

n 

ft 













o 

in 

o 

© 

o 

© 

© 

© 

o 

© 

o 

to 

o 

© 

o 




•* 

ft 

CM 

ro 

ro 

♦ 

<f 

in 

© 

© 

© 

ft 

© 


© 
eg 
o 
J © 
W o 
u o 
U « 

< o 

o 

tu o 

> • 

< o 


♦ 

N 

CM 

in 
<0 
«r 
o © 
< © 
O ft 

> • 
< n 
a © 

ft 


0 

1 

o 

© 


© 

© 

i 

o 

© 


© 

<y 



» 
































r> 





to 















o 














a 


a 

o 


a 



to 















eg 

o 

o 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a z 


a 

mm 


f* 



a 















• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

« 

• 

o o 


to 


D 



o 

o 

o 

O 

o 

o 

© 

o 

o 

o 

o 

o 

o 

o 


At 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 



X 

n 





• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

■ 

• 

• 

• 


















mm 

© 


U 



ft 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a 



w 

mm 

W 

w 

mm 

w 

w 

w 


w 


mm 

mm 

ft o 


J- 

ft 



o 































At ri 



• 







w 


w 


w 

w 

^m 


w 


mm 

w 


* 














o © 


a 

o 


ft 

o 




(Vi 

rvi 

■o 

© 

<» 

© 

< 

a 

© 

a 

u 

a 
















• • 


o 



t/> 



ft 


u. 

n 


ft 

H* 


1 

X 

1 

X 

1 

X 


ft 














© <y 


rr 






n 


7 

j 

<u 

m 

r 

r 


< 

_4 

< 

mt 

< 


J 
















a 






o 


►M 

tu 

X 

X 

a 

a 

a 

on 

o 

to 

o 

© 


UJ 














<g © 








a 


> 

> 

ft 

ft 



a 

* 

a 

X 

a 

* 


o 























• 

• 



































< 

o 


© 

< 

© 

<t 

© 


© 

< 

© 

«» 

© 


Af 

© 


© 

At 

© 

At 

© 

At 

© 

At 

© 

«T 

© 













N 

(V 

n 

n 

<T 

<T 

m 

to 

© 

o 





** 

(V 

n 

ro 

n 

At 

<t 

to 

© 

© 

© 




© 



ft 































2 Z 



© 



ft 


n 

n 



mm 


■mm 

mm 

mm 

mm 

mm 

mm 





mm 

mm 

mm 

mm 

A 


mm 

mm 

mm 

mm 

mm 

mm 


mm 

o o 



m* 



© 


o 


r> 













to 


© 















m4 



© 


© 


© 













o 


o 















ft 


































ft ft 

to 





© 


o 


O 













o 


o 












© © 

a 


© 

z 

a 

ft 


r> 


© 















© 












© © 

UJ 

*> 

© 

* 

o 

ft 


n 


<0 















© 












• • 

ft 


© 

o 

< 

ro 

a 

Oi 


w 












© 



n 












n © 

Ui 


en 

o 


►» 

tu 

N 


© 












a 

© 


mm 












o r> 

X 


• 

* 

z 

• 

ft 

(Vi 


O' 












UJ 

ro 

o 

© 

o 

o 

© 

a 

© 

© 

o 

O 

o 

© 

o 

mm © 

< 


© 


< 

10 

tu 

O' 


O' 













• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


a 



a 

i- 

© 

X 


o 

o 

o 

o 

o 

© 

o 

o 

o 

o 

o 

o 

o 

UJ 


o 

mm 

o 

o 

O 

o 

o 

o 

o 

o 

© 

O 

o 


< 



a 



< 

• 

• 

• 

• 

• 

• 

« 

• 

« 

• 

• 

• 

• 

• 

X 

i 


1 













a 



© 



c 

© 

o 

© 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

< 



w 

w 

mm 

w 

>— * 

•4 


mm 

mm 

w 


w 








< 

1 














ft 

M 


mH 












a »t 

z 



ft 



a 

w 


V 

w 



w 

mm 




«_, 

W 

_ 

«r 



II 












it a 

n 



z 




n 


mm 


m 

tn 

rv 

rQ 

< 

< 

n 

m 

u 

u 

a 

•rf 


z 












o o 




a 



ft 

X 

< 

u. 

o 

*- 

ft 

— « 

mrn 

1 

UJ 

i 

nr 

1 

UJ 


ft 


M 













© 



z 



7 

_■ 

A- 

z 

i 

UJ 

iii 

T 

T 

m 

T 


T 

m 

X 

t- 

-J 















-J 


© 

o 


Ift 

Ui 

tr 

-1 

•X 

UJ 

T 

T 

0 - 

a. 

c r 


a 

mrn 

a 

mm 

z 

•U 


• 












© 4t 

■> 


At 

a 


© 

a 

a 


> 

> 

ft 

ft 



a 

ft 

a 

ft 

a 

ft 

Ui 

o 


ft 












n cj 

a 


m 

r 


ft 

2 


• 

• 

• 

• 

• 

* 

• 





• 

• 

o 

• 

» 

* 

• 

* 

# 

• 

• 

■ 

• 

• 

• 

• 

♦ 

n <o 

o 


© 

o 

ft 

© 

UJ 

n 

© 

ro 

© 

n 

© 

n 

© 

n 

© 

ro 

© 

ro 

© 

z 

ro 

© 

ro 

© 

ro 

© 

ro 

<o 

m 

© 

ro 

© 

ro 

© 

• • 

a 

UI 

ro 

u 

z 

ro 

a 



m* 


eg 

(Vi 

ro 

r> 

•g 


to 

© 

© 

© 

UJ 



** 

•* 

eg 

eg 

to 

TO 

« 


© 

© © 

© 

© © 

a 

X 

© 

to 

< 

»0 

UJ 















a 















© At 


>— 

© 

-1 

ft 

© 















UJ 















ft 

u 

ft 

w 

to 

_J 

© 

z 

a 


mm 



mm 

mm 

mm 

<* 


4m 

mm 



o 



— * 


mm 

mm 

mm 

mm 

4m 



4m 




M 


© 

< 

a 

ft 


•M 

o 














to 















a 

a 

© 

X 

a 

o 


© 

o 

o 













o 















ft 

u 

• 


5 

• 


1 















1 















u 

a 

<r 


a 

mm 


o 

o 

o 













Q 















a 

a 

mrn 


a 

(VJ 


ft 

o 

a 













< 














z z 

j 


ft 



© 



o 

to 













o 














o o 

a 





mm 


ro 

o 

© 













® 






















ft 

o 

© 













in 






















CM 

o 

*4 













© 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

© 









o 














• 

• 




• 


• 

• 




• 


© 








ft 

o 

© 

o 

o 

o 

O 

o 

o 

o 

© 

O 

o 

o 


•4 

o 

o 

o 

o 

o 

© 

o 

o 

o 

o 

o 

o 

o 

O ft 









• 

« 

« 

• 

• 

• 

• 

• 

• 

• 

* 

• 

• 


t 














• • 



O 



n 


»0 


CM 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 





*4 

w 

mm 

mm 

mm 

*4 

s» 

w 

w 

mm 


O ft 



O 



ft 


1 















ft 














mm 



© 



* 



w 


w 

w 

Vr 

w 

w 

««* 

w 

mm 


w 

w 
















. © 


> 

o 


z 

© 


eg 

© 

ft 

Vt' 

OJ 

ft 

cu 

ft 

< 

< 

m 

© 

u 

V 


< 
















o 

o 


o 

mm 


Z 

© 

tu 

tr 

ft 

ft 



1 

a 

i 

ft 

/ 

t L 


ft 
















2 

o 


w 

ft 


*4 

<r 

-> 

< 

UJ 

L’J 

X 

I 

(V- 

2 

CVJ 

z 

eg 

Z 


J 
















a 

o 


© 

© 


•a 

X 


a 

X 

X 

a 

a 

a 

■" 

a 

MM 

n 



UJ 
















MM 

o 


-J 



a 

J 

> 

mm 

ft 

ft 



a 

> 

a 

> 

a 

> 


o 














ti. 


u 

10 


D 

© 




* 

• 













» 

• 

• 

• 

m 

• 

• 

• 

• 

• 

• 

• 


a 



© 


a 

eg 


(V 

ft 

(V 

ft 

(V 

ft 

CM 

ft 

rv 

ft 

cu 

ft 

Co 

ft 


eg 

ft 

(V 

ft 

(V 

ft 

W 

ft 

(W 

ft 

(Vi 

ft 

<g 

ft 

— o 


a 

• 


o 

• 




** 

w4 

eg 

tg 

n 

r> 

A? 


© 

© 

© 

© 





•M 

<y 

eg 

ro 

n 


At 

© 

© 

© 

© 

to 


a 

o 


a 

o 


v 

< N 





4m 


«M» 

MV 


PV 

mm 

p* 


mm 

A 

«*. 


«. 

4. 



p* 

PV 

mm 

-» 

p* 

p* 

PV 

p* 

D 

© 

ro 

mm 

At 

(V 











to 

© 

mrn 














o 

O 

o 

o 











© 

o 

O 













• 

1 

1 

1 












1 

1 













© 

© 

a 

o 

o 

o 











o 

o 

o 












a 

o 


© 

© 

o 











© 

w 

o 












z 

U) 

N 

H 

o 

o 











© 

© 

o 














(VJ 

© 

© 

ft 











© 

ro 

u> 












ft 


O 

ro 

At 

ro 











ry 

mrn 

ro 














n 

At 

© 

© 











© 

w 

o 

o 

o 

o 

o 

© 

o 

o 

© 

O 

o 

O 

o 


r> 

« 

© 

n 











• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

z 



w 

© 

♦ 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

ro 

At 

CJ 

o 

o 

o 

o 

o 

o 

© 

© 

© 

© 

o 


z 

• 

• 

* 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 















X 

o 


10 

eg 

© 

o 

o 

o 

o 

o 

© 

o 

o 

o 

o 

mm 

VP 

mm 

VP 

VP 

VP 

VP 

mm 

VP 

Vp 

VP 

mm 

VP 

mm 

u 



1 













X 

o 













ft 




*4 



VP 

*** 

vp 

VP 

*p 

v» 

VP 

— 

VP 


im 

a 












V* 

© 


r> ,J — ©< < Q £D u u 

IMUIUHI--" I I I 1 I I 

-ouiuuuii-f'i-n-n 
aou — xxaaaaciQcro 
a a < ft ft ft a a a a a a 

«* © v» © -• •Gw*'G~*>C’- t 'Ci*+'Q 
wHWNnnfltftfiinoo 


< Hi 
»- ft * 
j -j o 
ui u a 
o o 


mm © 


•A©*M©««A©**©p-©*pt0 

m *• n n n n < c tn m © © 


» © 

to © 

ft © 
© »o 
:> ♦ 
a 
x 


© 
© 
ft 
o 
a 
© 
© 
© 
a < 
ui to 
ft • 
o o 
a eg 


a n 
© 


© 

© 

o 

CM 

a 
© 
-J O 
< © 
m <f 

ft © 

M • 
2 © 
- © 
eg 
n 


(Vi 


J 

-J 

< 


ft 

ff 

0 

ft 

ift 

X 

ft 

z 

D 

0 
u 

1 

X 

u 

ft 

ft 

to 


< 

J 


70 


•*22 THRUST COMPUTE STEPS. 32 COAST COMPUTE STEPS 


>* w o o o r* h* o o ft © — ft 

<j»**»i»#**** 
ttooooor-oonoin<? 

kz © 0 0 ft ft ft ft 

< < 


cvo©©ooooooo©o 


ooc<?ooooooono 

n O - n 04 0 


(DiooriN^ nifloinoo 


D<f<f«ffV^OCOvDOOO 

<v<v<vrvcvcvcv 4444.400 


oooooooooooo© 


oooooooooooo© 
o o 


ootv^-onooooooo 


OO©© 0 CVOOOOOOO 

a 0 


Z 

o 


z 

o 


i~e:cvcv©©©ooo©''>© f ''o©©©cnffl© 0 rvn©o©©©©n©©©©©o 



G.*OOOOOOOOOOOOOOOO0000OOOOOOOOO©OOOOO 
2 D w (\J N N M N N NN N NWNf«^--NNN N N N N N N N N N ft| 

— 0 . 


© o o 

• • • 

o o o 


o 000 oo 44 * 4 ©©oo©©©©©©©©©©© 


000000000000000000000000 

MNMNNNNAINMAI 


z 

X 


z 

X 


z . 

44 

X • 


« C 1 O' n n in # 
I ***** * 
ft. h- f* N h P 1 * 
fl) <0 (0 CD iO * 


♦ (OwO’flo^fflNOocDtf e>Nojo*ao'Oooninio«o>* * ♦ * n o ^ in o id no o m o * ftNs&KNO^inento 

»•••** * * • « • • f ********* ■ • ft # » • t « • « * « ( t • 

« n 0 1 •* w ® * «N-'fflooN'Cin<f* nono>-h»«i^«-*^o>* * «- « ©h.©o©©©NNftnft*©©©oo©o 000 h- 

« r- CD o O 0 * Oiy'^COODCOON* K(DOCOD > (7'9 i OC'S)ft « « # OSOOOOffiOOOajft * ©©©0000©©©© 

• * • * « ft ft ft ft * 


in 
Ul 
J 

O ft 

Z *- 

<UJO© 00 ft» 
I CD O h N «» ft 


X 

< 

z 


X 

ft 

X 


X 

ft 

Z 


in © © © 


in 


D 

a 

x 


*ininft^ftftcvcvKmo. 4 tn«n©o*iniAnj©o 0 cvoo©**ftftNin©©® 44 ©<>n 0 »n**in©ft<*©o»>iocn©© 

#*«•••**•••»••••»* •••••»*(»«ft«ft» ••t»»**i*«*ftft«*tt***i*** 

ftioo | ^^<D*n > N ,, ) 4 <t(\jo'noo»noN<o<nn(v 9 i>*oiftftftftiDNOAiK<\iNHO(\ioftftcooH^<C'i)>*Ho 09 

♦ J CDf 0 ino* 00 O®‘D©NNN*t*N 0 ON- — nNDHiO»ftft«<lO«l/)in' 0 ' 0 'OO^r)»«OSS«<r)DMOI 9 N <0 

* | 44 4 * .4 # * # f | | *4 *4 r 4 44 

*11* » III *•** *» III 


z 

*-e 

z 


ft 

z 


u 

z 

D 

IL 


o 

Ui 


M 

ft 

* 

* 

CV 

ft 

OJ 

© 

o 

ft 

© 

44 

N 

•4 

44 

o 

ft. 

© 

© 

44 

* 

CD 


© 

o 

© 


© 

© 

o 

(VI 

* 

* 

« 

’* 

n 


© 

© 

4 * 

0 

© 

© 


© 

© 

ft 

ft 

© 


n 

* 


ft 

ft 

Cd 

•4 

0 0 

• 

• 

ft 

#• 

• 

• 

• 

• 


« 

• 

• 

• 

• 

• 

• 

• 

• 


• 

ft 


• 


• 


• 


• 

• 


* 

« 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 


• 

• 

ft 

ft 


• 

• 


• 

• 

• 

• 

• 

• • 

f*- 

O 

ft 

ft 

m 

in 

N 

N 

ft 

ft 

© 

© 

<1 

© 

© 

© 

s 

© 

if) 

© 

ft 

© 

n 

© 

ft 

r* 

m 

© 

N 

© 

o 

* 

ft 

ft 

* 

© 

N 

N 

N 

N 

© 

© 

© 

© 

eg 

44 

* 

ft 

O' 

O 

O 

© 

O' 

O' 

O' 

© 

© 

© © 

N 

<£> 

ft 

ft 


© 

© 

© 

© 

ft 

© 

© 

N 

K 

r* 

N 

ft- 

ft* 

ft. 

© 

ft 

N 

© 

© 

© 

© 

© 

© 

© 

© 

© 

* 

* 

* 

* 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

ft 

ft 

N 

© 

© 

© 

© 

© 

© 

© 

© 

© © 

1 

1 

ft 

* 






• 

1 

1 

1 

1 

• 

1 

1 

( 

1 

1 

ft 











ft 

ft 

ft 

ft 

1 

1 

1 

1 

1 

1 

I 

l 

I 

l 

1 

ft 

ft 











z 







X 






Z 














X 










Z 

















X 











< 






•4 














< 










tat 

















< 




X 







z 






z 














z 










z 

















z 





u 

UJ 

-J 

u 

v> 


r o o o o o 
u - t 

>- >- o a o o o i 


w “> • 
Ik. ft 


i *-• 

© 

o 

o 

o 

ft 


© 


o 

O 

o 

o 

o 

o 

o 

© 

o 

ft 

O 

o 

o 

O 

o 

o 

o 

O 

© 

o 

o 

o 

4 * 


o 

o 

o 

o 

o 

o 

o 

o 

44 

© 

O 

44 

© 

44 


O 

o 

o 

o 

o 

o o 

44 

» o 

-* 

o 

o 

o 

•J* 

o 


o 

o 

o 

o 

o 

o 

o 

o 

Y 

o 

44 

O 

o 

o 

O 

o 

o 

O 

O 

— 

o 

o 

o 

o 


o 

o 

o 

o 

o 

o 

o 

o 

o 

44 

O 

o 


o 

o 

o 

o 

o 

o 

o 

o o 

o 

1 o 

o 

o 

o 

Ct 

O 

o 

a 

o 

a 

o 

O 

o 

a 

o 

o 

O 

o 

o 

o 

o 

O 

o 

O 

a 

o 

o 

o 

o 

o 

o 

o 


o 

o 

o 

o 

a 

o 

c> 

© 

c> 

o 

o 

o 

o 

a 

o 

o 

o 

o 

o 

o 

o o 

o 

. m 

© 

r> 

© 

o 

in 

<n 

rh 

o 

<» 

K 

o 

r» 

m 

*♦ 

© 

ru 

ft 

cv 

cv 

•4 

K 

0 

© 

M 

© 

© 

o 

44 

cv 

ft 

rd 


o 

ft 

© 

© 

44 

© 

44 

© 

0 

* 4 , 

o 

44 

cv 

K 

cv 

0 

ri 

ft 

ft 


4 . © 

© 

Al 

© 

© 

o 

o 


o 

m 

lA 


4 * 

ft 

m 

+ 

© 

*n 

OJ 

til 

ft 

© 


nr 

N 


»n 

rv 

in 

O 

in 

K 

o 

0 

r> 

— 

o 

Al 

»A 

N 

© 

** 

-4 

cv 

© 

o 

«n 

Al 

«n 

nt 

cv 

m 

m 

m 

*- 

K- © 

m 

© 

•p 

CM 

tu 

Cd 

l 

k 

til 

© 

- 

m 

© 

n 

n 

til 

cv 

cv 

1 

1 

l 

- 

- 

CV 

CV 

CV 

CV 

cv 

- 

CV 

1 

I 

1 

1 

CV 

1 

o 

- 

cv 

cv 

cv 

IV 

cv 

rv 

(V 

ft 

« 

cv 

1 

N 

1 

CV 

1 

- 

cv 

cv 

c\i 

cv 

cv 

- 

44 44 

0 


7 


X 


U. 

Z 

Z 





X 




U- 

z 

z 




X 





b. 


z 



Z 





X 





b. 

z 


7 





X 




z 


O 


< 


U. 

*• 

o 





< 




u. 

*4 

o 




< 





IL 





u 





< 





u. 

44 


U 





ft 




u 




X 


u 

z 






z 




o 

z 





z 





o 


z 








% 





O 

X 







z 







© 

© 

N 

© 


CM 

IV 

tv 

ft 

ft 

N 

o 

0 

0 

© 

o 

ft 

N 

N 

K 

K 

N 


K 

0 

ft 

K 

© 

0 

o 

o 

© 

ft 

cv 

© 

n 


0 

© 

© 

K 

0 

m 

ft 

ft 

ft 

til 

o 

0 

© 

rv 

© © 

© 

© 

n 

ft 

© 

© 

K 

© 

© 

© 

® 

r* 

►- 

K 

r* 

© 

© 

N. 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

0 

0 

0 

0 

o 

o 

0 

0 

0 

© 

© 

© 

K- 


® 

© 

0 

o 

o 

o 

o 

0 

0 

m 

© 

© 

© © 

0 

O 

o 

O 

o 

O 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

O 

O 

O 

o 

O 

o 

o 

- 

44 

o 

o 

o 

o 

o 

o 

O 

o 

O 

O 

o 


- 

44 

44 

o 

o 

o 

o 

o 

o o 

o 









X 






z 







X 


z 








X 









z 







X 





Z 











< 






M 







< 


M 








< 









•4 







ft 





44 











z 






z 







z 


z 








z 









z 







Z 





z 



© 

© 

til 

O 

© 

n 

n 

© 

in 

N 

ft 

© 

0 

- 

© 

© 

(V 

© 

*4 

CV 

cv 

n 

rv 

© 

© 

44 

© 

0 

n 

rv 


cy 

n 

- 

© 

© 

ft 

© 

ft 

s 

© 

© 

© 

© 

44 

© 

0 

© 

4 

0 

r> 

© 

r* 

K f* 

44 


J, 

ft 

© 

h* 


n 

ft 

© 

© 

o 

o 

o 


o 

0 

© 

cv 


o 

o 

o 




CV 

IV 

44 




44 

*4 

cv 

© 

ft 

© 

© 

K 

N 

c* 

ft 

n 

0 

© 

r* 

N 


cv 

CV 

© 

© 

© 

© © 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

N 

H 

** 

N 


© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© © 

© 

z 













X 







z 






X 




z 









X 







z 







■ X 















< 













< 













< 







*4 







ft 


z 













z 







z 






z 




z 









z 







z 







z 


© 


o 

o 

© 

o 

tu 

n 

ft 

O 

0 

0 

© 

© 

© 

o 

© 

© 

K 

© 

ft 

© 

cv 

© 

0 

ft 

o 

rv 


Cd 

© 

0 

o 

n 

© 

© 

ft 

n 


0 

K 

o 

0 

© 

© 

© 

© 

© 

o 

o 

o 

o 


44 f 4 

*4 


© 

m 


ft 

n 

© 

N 

© 

o 

4 * 


*4 

•4 

o 

o 

r» 

© 

© 

© 

© 

© 

© 

© 

© 

N 

© 

0 

o 

o 

o 

o 

44 

44 

44 

44 

44 

•4 

44 

o 

o 

o 

0 

0 

0 

0 

0 

0 

o 

o 

o 

o 

o 

o o 

o 

© 

© 

© 

© 

© 

o 

o 

m 

o 

O 

o 

o 

o 

o 

o 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

0 

0 

0 

0 

0 

0 

o 

o 

o 

o 

o 

o o 



< «i N 


a 

o. 


Z X 

o < 

z 


u. z 

IL - 

o z 


z 

D 


o © © ft © 


z . 
o 


i- in 
< *1 
u o 


z 

D 

Z I 
Z . 
O I 

u : 


-J : 
o 
in i 


l/)«*M*4*4OOOOO0OOOOO' 

o 


-X 

ft 


0000000000*4 

z * 


I O O O O O O I 


o 

o 

© 

K 

IV 

© 

© 

© 

0 

0 

rv 

© 

© 

© 

© 

0 


cv 

r> 

© 


© 

0 

© 

0 

- 

f* 

0 

© 

o 

cv 

44 

O 

© 

© 

44 

© 

0 

o 

ft 

cv 

0 

© 

© 

(V 

© 

0 

© 

0 

© 

N 

r> 

© 

© 

© 

© 

*4 

•ft 

0 

© 0 

P 4 



K 

N 


© 

© 


© 

© 

ft 

0 

© 

o 

0 

ft 

m 

© 

cv 

44 

o 

0 

« 

© 

© 

n 

0 

© 

ft 

J 

0 

o 

ft 

© 

o 

© 

© 

© 

n 

0 

© 

44 

o 

cv 

o 

0 

© 

© 

44 

© 

© 

o 

CV 

n 

© 

O 

n 

ft 

ft 

© 

© 

0 








© 

© 

0 

© 

N 

0 

0 

44 

© 

© 

© 

N 

0 

0 

n 

w> 

N 

o 

O 

44 

ft 

© 

© 


© 


r> 

n 


© 

© 

h- 

0 

o 

44 

rv 

m 

© 

ft 

K 

© 

44 

n 

© 

© 

© 

0 

o 

o 

O 

W 

cv cv 









*4 


4 * 


cv 

(V 

IV 

cv 

cn 

r> 

r> 

n 

m 

n 

n> 

ft 

ft 

ft 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

N 

f* 



r* 

K 

r* 

N 

© 

© 

© 

© 

© 

© 

0 

0 

0 

0 

0 0 

0 


© 

0 

0 

0 

© 

K 


ft 

© 

© 

f* 

0 

© 

0 

44 

© 

0 

© 

rv 

© 

ft 

N 

cv 

© 

rv 

o 

© 

ft 

cv 

o 

h* 

o 

CV 

o 

cv 

© 

cv 

© 

ft 

ft 

© 

© 

© 

© 

ft 

© 

© 

4 * 

© 

© 

K 

44 

© 

ft 

© 

© 

rv 

N 

4 * 

ft 

o 



N 


n 

© 

K 



© 

© 

ft 

ft 

*4 

© 

f* 

© 

c 

m 

© 

f* 

ft 

O 

n 

0 

© 

n 

ft 

© 

© 

cv 

ft 

K 

0 



ft 

0 

© 

0 

© 

C* 

cv 

N 

© 

© 


© 

<v 

© 

ft 

© 

N 

o 

CV 

n 

© 

© 


© 

O 

o 

o 

o 

© 

ft 

© 

0 

© 

© 

0 

© 

© 

o 

cv 

ft 

— 

© 


0 

Y 

© 

a. 

© 

o 

ft 

0 


o 

0 

o 

© 

© 

© 

— 

— 


o 

cv 

© 

N 

© 

0 

ft 

r* 

ft 

0 

n 

© 

© 

© 

CV 

0 

0 

•* 

o 

© 

*"■ 

ft 

Cd 

© 

© 

4 * 

© 



© 

0 





© 

0 

© 

o 

•ft 

© 

n- 

CV 

© 

© 

© 

© 

© 

© 

J, 

o 

0 


f* 

© 

ro 

0 

© 

0 


© 

© 

n 


m 


n 

- 

© 

© 

© 

rv 

ft 

44 

n 

0 

c 

© 

© 

0 

o 

K 

n 

© 

© 

© 

c 

® 

cv 










© 





© 

r» 

44 

44 - 

cv 

n 

© 

© 

o 

cv 

m 

© 

© 

CV 


o 

44 

rv 






4 , 

44 

n 

© 

© 

© 

© 

0 

o 

-4 

ft 

<3 


0 

O 


44 

© 





0 

0 











44 


(V 

cv 

CV 

rv 

rv 

r> 

n 

© 

n 

n 

m 

ft 

ft 

ft 

ft 

ft 

ft 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

t- 

N 

n 

t •* 


N 

© 

© 

© 

© 

tt 

© 

© 

© 

© 

© 

0 

o 

ft 

© 

ft 

ft 

© 

»* 

© 

ft 

ft 

© 

ft 

© 

© 

« 

ft 

© 

ft 

© 

« 

ft 

n 

© 

© 

© 

ft 

ft 

© 

ft 

« 

© 

ft 

ft 

ft 

© 

© 

© 

© 

© 

© 

© 

ft 

ft 

© 


ft 

© 

ft 

ft 

© 

© 

© 

© 

ft 

ft 

© 

ft 

© 

© 

© 

ft 

© 


_J u 
U z 
tu o 

-I 


tu 

z 


Ul 

2 


nff 


?ooft 


a 1 


71 


MISSI3N SCHEDULE 


© © 

• © © 

« — — 

z • • 

o O' c 

mi < « 


© DO 

3 <* < 

*• N M 

O O CO 

< o o 

0: <f < 

© o 

o o 


0 © 

1 I 
o o 

*- m r** 

o — • <NJ 

O -4 O' 

>• o 

CM O 
© (0 


o o 

L- Nm 

n «t © 

o mo 

x n 


0 o 

1 1 

Q O 
© © 
o o 
> o>o> 


n 

© 

© 


0 

> © 

© 

► 



<r 0 

« 

K 

< 

z 


< ♦ 

a 

m* 

. 0 

Ul 


r - 

K 



wa 


u 



© 

u ~ 





0 

N- u 

(ft © 




© 

tL UJ 

© © 



c 0 

1 0 

u. (/) 

< N 



n 

« 

UJ V 

z < 

0 

© 

0 

w 

0 z © 

• 

u 0 

V 

CM 

•4 

U X • 

»- to 

— © 

< 

• 

© 

— 0 

to © 

a n 

O N 

© 


V 

Z h* 

a © 

w ffi 

© 


O D 


to • 

© 

* 

u 

z 


0 

UJ — 

r> 


UJ 



z 0 


to 

M 



K « 


z 

u 



K © 





h- 


tL 




a 0 


O N 


© 



O 



a 0 




w 0 


a 


N 



O 

»- 


UJ 0 




to 0 


0 


Z 



© 

r 


0 


UJ 


© • 


UJ 


w 



© 

© 


to to 


a 


W O 


Q 


0 



• 

M 


•0 


3 


0 




a • 



© 

J 


• 


K O 


a 0 




Z 0 


© 

to 

IL 


© 


U • 


© n 




■4 


Ul 

M 





3 O 






> 


r> 

N 





a 









vs 





K 




© 


a 



•4 





© 

v. 



> 


Of 







«■» 


a 


© 

< O 


< 




V- O 



O 


< 

■Ml 

Z 

O CO 






Z K 



¥ 


z 

<u 

O 

w CM 






0 to 





z 

« 

►4 

© 




to 

II 

© 

© 




3 

• 

H- 

UJ c> 




© 


© • 

a 


V 


© 

U <t 

M 

z • 




© 

© 

a n 

UJ 


a 



to O 

O 

— 

© 



O' 

K 

3 © 

K 


< 


to 

© 1 

z 

H — 

Z 



4 

Z ' 

O © 

tu 


z 

Ul © 

u 

V O 

0 


a 



* 

(0 

Z * 

z 


z 

© K 

z 

Z 0 

0 

z 

M 

* 



*■4 

« 

< 


3 

< m 

< 

w CO 


ff 

L- 

« 



u 

O 

a 

a 

© 

¥ S 

z 

< 

Ui 

D 




n 

M 

<D 

< 

0 


Z • 

or 

— © 

u 

© 

O 

UJ ^ 

V 


tL 

to 

a 

L* 

© 

< to 

0 

3 © 

z 


z 

© M 

a 

z 

IL 

<0 UJ 


u 

© 

K O 

u. 

< © 

< 


a 

V. CvJ 

< 

0 

to 

• 

z 

< 

< 


a 

• 

z 


u 

IM © 

z 

■M 

O 

© < 

to 

IL 

z 


UJ 

-4 «\f 

a 



* O 

z 

t- 

u 

0 

L- 




f L 


0 

• ^ 

_j 

* O 

3 

< 



© 

10 

z 



O 

a 

Z - 

< 

z • 

© 

z 


• z 

>* 

a 0 

UJ 


z 

O 

or 

'w <t 

> 

¥ © 


M 


«- < 

© 

3 • 

►* 


UJ 

< 

UJ 

m* 


w CM 

UJ 

-J 


CD — 


h- 0 

© 


K 


a 

H © 

a 


0 

u 


© _J 

Z 

u 

► 


© 



© © 

a 

to 

z 

. z 


~3 

O 

0 

© 


> 


D 

3 - 

< 

U 

< 

*4 


"I 

M 

rr 



© 


z 

a 



z 



• 

© 

L* 

z 

♦“ K 


Mi 

< 

T <T> 

0 


fT 

-J 


0 

J 

© 

0 

Z m 

z 

z 


L* • 

z 


0 

< 


CM 

© 



< N 

a 

W m* 

> 

O 

< 


u. 

z 



a 


© 

-J 0 

M 

© 

a 

X 



a 

M 


»- 

0 


-1 

-1 • 

© 

— © 

0 

<r 

UJ 


to 

IL 


y 

a 


3 

Ui -« 

-1 

3 - 

►- 

Z 

a 

*+ 

0. 



0 

a 


a 

a n 

D 

< N 

u 


3 

O 






a 

0 

O to 

a 

to 

UJ 


L- 

© 




11 

u 

0 

a 

a *+ 

0 

w* • 

3 


a 

© O' 






*- 

a 

CL 

a 

w 0 

< 


< 

V «t 




0 

a 

u 



a 

cc 

CC 


a 

Z to 




< 

H 

< 0 

0 



•T 

b- 

M 

UJ 

0 CVJ 





u 

IL © 

*4 


© 



* 

0 

© 





UJ 

r> 

a 




UJ 

¥ 


a • 





• -J 

to 0 

►* 


a 


z 

- a 


z 0* 





<C O 

- a 


u 



a 

a 0 

> O 

z 


►- 0 

¥ 

mi 

fc 0 

UJ 



V- 

UJ o* 

© 

0 


Z N 

Z 

UJ 

Z 0 

mi 

ft 


X 

ft «* 

a 


or 

0 n 

< 


< © 

UJ 

¥ 

PJ 

UJ 

0 © 

UJ 

© 

3 

0 

H 


J w 



CM 


a • 

0 

© 

< 

© • 



a • 



n» 


0 


IM 

►- 

a © 



0 


»*» 

• 


X CM 


z 

Z 

3 ® 



a w 


O 

0 


< 



UJ 

O O' 



10 © 


a 

CM 


z 


u 

O 

X n 



ft 


< 





** 

\ 




0 


H 





L- 

UJ 

0 

V 


a 


w 




M 

a 


to 

¥ 




a 




O 

*4 

M 

© 

S 








UJ 

-j 

•4 

® UJ 

O O 








0 0 

0 

M 

• L* 

¥ O 








44 © 

UJ 


<VJ < 

w O 





* 


3 

O 


z 

-4 D 

O 




mm 



< N 

O O 

a 

< 


h • 


j 0 


ft 

O 


w © 

Z • 

0 

*- 

•Z 

m 


< © 


¥ 

O 


0 

*- © 



< 

< 44 




w 

© 


H © 

© 

f- 

H 

N - 

z 


*- © 



© 


© 0 

¥ 

3 


O' mi 

a 


*4 • 



• 


•- CM 

a 

O 

© 

~ 3 

j 


z © 


a 

O 


O © 

< 


*4 

3 

< 


- © 


¥ 



• 

a 



• 



ti 


© 



z 0 



UJ 

*4 



n 


Z 



•4 



A 

N 








z 



0 






a 






. 44 

a 











r 

a 

M 










UJ 

< 

ft 








M 


> 


¥ 








O 



p 

S 




M 



3 

UJ 


Z 


3 0 




ft 

O 


« a 

0 K 


u 

0 

¥ O 




¥ 

O 


w O* 

w n 


72 


B. COMET RENDEZVOUS MISSION 

The objective of this mission is to deliver maximum payload to the comet 
Tempel II, rendezvousing at perihelion of the comet's path in the 1988 apparition. 
For the case shown, the arrival date is fixed at September 16, 1988; the launch 
date and launch excess speed are optimized to yield maximum payload for the 
nominally 4-year class of solutions. The launch vehicle assumed is the Titan 
III E/Centaur. The electric propulsion system parameters are representative 
of projections of the SERT III spacecraft.. Specifically, the reference power is 
8.671 kw, the specific impulse is 2900 seconds, the efficiency is 0.63376 
and the specific propulsion system mass (arrays plus power conditioning and 
thruster subsystem) is 30.285 kg/kw. 

This example illustrates the use of the solar array degradation option. 

This option is invoked by setting the characteristic degradation time TPOWER 
to a value less than 10 10 . In this case, it is set to 7121 days which means 
that if the arrays were situated at 1 AU and oriented normal to the sun line for 
this amount of time, the power developed by the array would degrade to l/e of 
its initial power output. The input MPOW = 1 forces the arrays to be oriented 
normal to the sun line throughout the mission. Setting the triggers X30(2) and 
Y30(2) to 1 results in the optimal adjustment of the trajectory to accomodate 
the degradation. The complete set of inputs for this case follows and the resultant 
program output begins on the next page. 


&MINPUT X1(2)=1.DO,X2(2)=1.DO,X3(2)=1.DO,X4(2)=1.DO,X5(2)=1.DO 
X6(2)=1.D0,X11(2)=1.D0,X13(2)=1.D0,X15(2)=1.D0,X7=1.D0 
Y1(2)=1.D0 / Y2(2)=l.n0 / Y3(2)=l.n0 / Y4(2)=1.D0,Y5(2)=1.D0 / Y6(2 ) =1. DO 
Y13( 2 ) =1. DO, Y15( 2 ) =1 .DO, ALPHAT=15 . 285D 0,MBOOST = 1 5, MTMAS S = 3 
MOPT2=3,MOPT3=43,MYEAR=1988,MONTH=9,MOAY=16,HOUR=17. 22D0 
X12 = 2 . 84 3928 5D4 , B I =. 6 3376D0,D I =0. DO, CTANK = . 1D0, Y 11=8. 6 7 IDO, 3. DO, l.D-3. 
TP0WER = 712 1. D0,M P0W = 1, X30(2)=1,D0,Y30(2)=1.D0 
Xl= 5. 072998588224D 00, X2= 3 . 3501 835 72 126D 00, X3=- 2. 1 59 80814 2 3 89D 00 
X4=-2. 069314933312D 00, X5= 4 . 56 524 774 1305D- 00, X6= 3. 317322088503D-01 
Xll = 1.647979654151D-04,X13= 6 . 718039 798493D 03, X15=-l. 5176252844 70D 03 
X30 =-4.6574062226460- 03 SEND 
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NOMINAL TRAJECTORY 1 (TOTAL I) INHIBITOR IS 5.82080-11 
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BOOT TRANSFER ANGLE BETWEEN EARTH ANO TEMPEL 11(19881 IS 10.9153 DEGREES 


CASE t (CONVERGED) PERFORMANCE S JMMARY 
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c. MULTIPLE BALLISTIC SWINGBY MISSION 

This computer run demonstrates how the program may be effectively used 
in the investigation of all-ballistic missions (with electric propulsion absent)., 

In general, trajectories are forced to be all-ballistic by setting the pro- 
pulsion-time adjoint variable X8 to a "large" negative number (X8 = - 1.D3) 
with respect to the mass ratio adjoint variable X7 = 1.D0. These settings (to- 
gether with the reference thrust acceleration Xll = l.D-4, jet exhaust speed 
X12 = 1.D0, and non-zero primer vector components X2 = .IDO at the launch 
planet and X42 = „1D0 at the first intermediate -target) guarantee that the 
electric propulsion thrust switch-function will be maintained negative, yielding 
ballistic flight at all times. The electric propulsion specific masses ALPHAA 
and ALPHAT are also set to zero, so that the net spacecraft mass equals the 
launch vehicle payload. 

The particular ballistic mission simulation demonstrated here involves 
a March 6, 1985 launch of a 1,635 kg payload by an Atlas/Centaur launch 
vehicle (using the launch vehicle selection default -value of MBOOST = 0; actually, 
changing MBOOST alters the payload mass computation but not the C or the 

O 

trajectory). 

The primary target is the Earth, specified by MOPT3 = 3, and there is 
one intermediate target, the comet Giacobini-Z inner, specified by MOPTX = 41. 
The launch occurs fifteen days before (X15 = - 15. DO) the reference date 
(MYEAR = 1985, MONTH = 3, MDAY = 21), and the spacecraft passes Giacobini- 
Zinner 174.22 days (X48) after the reference date and arrives back at Earth 
353 days (X16) after the reference date. 

Ten days after passing Giacobini-Z inner a deep-space burn (of 188.9 
m/sec) makes possible the re-targeting to Earth. This is accomplished by 
inputting TDV = 1. 00010D5 together with X64(2), X65(2), and X66(2) = 1.D0. 
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The deep space burn velocity-increment (X64, X65, X66) initial guess is the 
zero vector, by default. The deep-space burn could just have well occurred 
before arriving at Giacobini-Zinner (e.g., TDV = 2.00005D5). The iteration 
sequence consists of a 6x6 hunt, with the independent variables being the initial 
heliocentric velocity at Earth (X18, X19, X20) and the deep-space burn velocity 
increment (X64, X65, X66), and the dependent variables being the position target^ 
ing at Giacobini-Zinner (Y41, Y42, Y43) and at Earth (Yl, Y2, Y3). The setting 
MAXHAM = 0 is made to avoid an unwarranted BAD HAMILTONIAN warning 
message due to the presence of the deep-space burn. 

The detailed printout of the iteration sequence is omitted by setting 
NPRINT = 3 (compared to the default value NPRINT = 7), and three extra lines 
are added to each print -block by setting MPRINT = 2. 

Finally, the multiple ballistic swingby option is invoked by means of 
MOPT4 = -3, 42, which directs the spacecraft to swingby the primary target 
(Earth) in such a manner as to re-target to Earth (MOPT4(l) = -3), and then, 
at the second Earth passage, to again swingby in such a manner as to target to 
the comet Borrelly (MOPT4(2) = 42). Furthermore, the first Earth swingby is 
specified as unpowered (MSWING(l) = -1) having T2(l) = 480 days as the initial 
guess of the Earth-to-Earth transfer time, and the sec raid Earth swingby is 
specified to be powered (MSWING(2) = -5) with a specified transfer time from 
Earth to Borrelly (T2(2) = 127 days) and using the input Earth-departure helio- 
centric velocity initial-guess "XSWING(1, 2) = . 62D0, . 91D0, 0.D0. Both swing- 
by iterations converged, and the powered-swingby incremental speed turned out 
to be -13.7 m/sec, the minus sign denoting a braking burn. The default value 
of TGO =•- 1.D0 resulted in the printout of several pages describing the tra- 
jectory segments following the time of primary target passage. 

It should be noted that the trajectory simulation would have terminated 
at the primary target (first Earth encounter) if MOPT4(l) were zero, and would 
have terminated at the second Earth encounter had MOPT4(2) been zero. The 
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mission, as simulated, thus consists initially of a eomet-flyby, followed by a 
second comet-flyby making use of a double Earth-swingby to reach the second 
comet. 

The inputs for this computer run are listed below, and the resulting pro- 
gram output is displayed on the following pages. 


&MI NPUT X7 a l.D0,X2 B .lD0,X42 a .in0,X8=-l.D3, ALPHAA = 0 .D 0, ALPHAT = 0. D 0 

X6 4(2)=‘l.n0 / X65(2)=l.D0 / X66(2>=l.D0 / X18(2)=l.D0 / Xin(2) : O,D0 < X20(2)=l.D0 

Yl(2)=l.D0 / Y2(2)=l.n0,Y3(2)=l.n0 / Y41(2)=l.n0 / Y!t2(2)=l.Dfi / Yl*3(2)=l.n0 

MYEAR=1985,MONTH=3,MOAY=21 / MOPT2=3 / MOPT3=3 / MOPTX=41 / Xll=l.D-4,X12 a l.D0 

TDV=1.6ooiOD5,X15=-15.DG,X16=353.DO,X48 a 174.22DO 

X18 — 8. 32760-2, X19 — 1. 0034D 0, X2 0=-2 . 01430-3 

T2 a 4 80 . DO, 12 7. DO, XSW { NG( 1, 2) =. 6 200, . 91D0, 0 , DO 

MO PT4 =- 3, 4 2/MSW I NG=- 1, -5,MPR I NT=2,MAXHAM =0, NPR I HT=3 ft END 
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